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I. INTRODUCTION 

The G. E. Image 100 System provides the framework around which to 
build a versatile tool for the development, testing, and evaluation of 
image analysis techniques used in support of current and future Earth 
Resources Programs. In order to effectively support current applica- 
tions, the existing Image 100 System requires augmentation through the 
acquisition of additional hardware, the design and implementation of 
application-oriented software, and the development of procedural 
guidelines for efficient System utilization and management. 

An augmentation plan for the Image 100 System should be developed 
to provide the following capabilities: 

• A batch-processing mode for testing the effectiveness 
and applicability of newly-developed data analysis 
algorithms and techniques. 

« A means whereby interactive image analysis techniques, 
which are directly applicable to Earth Resources data- 
processing requirements, can be developed, 
s A means whereby data-analysts and analyst-interpreters 
can effectively interact. 

Developmental work aimed at providing the above capabilities can 
be divided into the following general areas and associated tasks: 
Procedures Manual Development 
Tasks 

1. Procedures Manual Outline 

2. Development of the Procedures Manual by completing 
each section of the outline 



Operating System Development (RSX 11D) 


Tasks 

1. Implement Operating System which provides batch/ 
conversational and interactive modes 

2. Maintain Operating System 

3. Modify Operating System to accomodate advanced 
file-handling techniques 

Image Analysis Library Development 
Tasks 

1. Identify computer programs which should he part 
of the Image Analysis Library 

2. Document, prioritize implementation of, and 
procure programs for Image Analysis Library 

3. Design Library programs so that they can be 
used effectively in interactive mode 

4. Implementation of Image Analysis Library 

5. Maintain Library 
Interactive System Development 
Taske 

1. Augment Interactive System which operates under 
Operating System (RSX 11D) with minimal modification 
of Operating System 

2. Augment advanced file-handling capability for both 
batch/conversational and interactive modes 

3. Interface with personnel developing Operating System 
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4. Maintain/Modify Interactive System 

5. Study hardware/software configuration for advanced 
operational system 

Interactive Capabilities Development 
Tasks 

1. Determine those interactive capabilities needed 
to support various Earth Resources Projects 
(LACIE, FAP, etc.) 

2. Implement software to provide interactive capabilities 
under Interactive System using Image Analysis Library 


II. SUMMARY OF WORK PERFORMED 


The work performed under the subject contracts is a portion of 
that described in the Introduction — Procedures Manual Development and 
Image Analysis Library Development , namely; 

® Preparation of an Image 100 Procedures Manual Outline 

® Identify and prioritize implementation of computer programs essential 
to earth resources data processing requirements (FAP, 

RAP, LACIE, RT&E, etc.) 

The outline of an Image 100 Procedures Manual was developed which 
sets forth guidelines, of both a general and specific nature, that 
provide a basis for the preparation and updating of an Image 100 
Procedures Manual, The scope of the outline was limited to definition 
of general features of a procedures manual together with special 
features of an interactive system. The outline is organized in a 
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"task to be performed" fashion, and the topical order is intended to 
be a logically ordered "table of contents." This work is discussed 
in more detail in the next section and Appendix A. 

Computer programs were identified (source, purpose, 
availability, application, etc.) which should be implemented as part 
of an applications-oriented library for the Image 100. In addition, 
each Identified program was assigned an implementation 
priority. Although the Image 100 System is an interactive system that 
will also operate in the batch mode, identification of computer 
programs was based only on the requirement that they be utilized in 
the batch processing mode. Computer program identification was 
limited to those having direct or anticipated application to earth 
resources data processing requirements. It should be noted that 
many of these programs are essential to the development of inter- 
active data processing capabilities not currently available on 
the Image 100 System. The Image Analysis Library is discussed in 
more detail in Section IV and Appendix B. 

III. PROCEDURES MANUAL DEVELOPMENT 


This portion of work is concerned with developing an outline of 
a Procedures Manual for the Image 100 System. By Image 100 System 
we mean the PDP 11/45 general purpose computer, special image 
analysis hardware (the G.E. Image 100) and all associated software. 
Even though the scope of this task is mainly limited to the develop- 
ment of a Procedures Manual outline, we did author those sections 


dealing with updating standards, program documentation standards, and 
library documentation standards. Also, we authored a JCL Cookbook 
for FORTRAN IV users. 

The purpose of a procedures manual is to provide information 
necessary to effectively use the system. This includes descriptions, 
operating procedures, maintenance procedures and so on (much of 
which is furnished by the manufacturers of the hardware), but must 
also include standards for documentation, maintenance and augmentation 
of system hardware and software. We feel that work performed in 
furnishing suitable documentation for a given program is more than 
offset by the long term reduction in effort in use of the program. 

Thus we pay particular attention in this outline to setting a 
flexible but high standard for documentation for programs at each 
developmental level. 

In formulating documentation standards, an effort was made to 
utilize existing widely-accepted documentation guidelines. While not 
properly a part of a procedures manual, the information we located 
may be of use to some programmers in their attempt to understand and 
follow the standards for documentation set here. We found the 
following documents to be particularly helpful; 

NASA Handbook, NHB 2411.1 July 1971 

"Computer Program Documentation Guideline" 

Goddard Space Flight Center, X-540-7 2-114 Feb. 1972 

"A Programmer's Guide to the Goddard Space Flight Center 

Computer Program Library" 

DOD Manual 4120. 17-M Dec. 1972 

"Automated Data System Documentation Standards Manual? 


Some sort of committee or working group should be formed to decide 
what is acceptable (well-documented, tame) enough to be implemented into 
the Image Analysis Library or into the Interactive system. At the 
present time, the Image Analysis Library Is a collection of forms which 
identify programs of possible use and a list of functional requirements. 
As a useful program is located and implemented, it can pass through the 
following stages : 

1. Identified (by a Computer Program Identification Form) 

2. Compatible Library Unit (card deck and documentation of 
a program compatible with the PDP 11/45 system) 

3. Specially-designed Library Package (a collection of 
compatible units assembled into an efficient package) 

4. Implemented Library Program (a completely documented 
program of general interest with object code kept 
on the system disk) 

5. Interactive Unit 

Library management information in the Goddard manual mentioned 
above may be a little elaborate for the Image Analysis Library, but it 
has not been ignored. However, we feel it would be premature to set 
final management standards for a library which does not yet exist. 

IV. IMAGE ANALYSIS LIBRARY DEVELOPMENT 

This portion of work is concerned with the identification of, 
and establishing implementation priority for various software units 
which will comprise an Image Analysis Library to support Earth 
Resources data processing requirements. By Image Analysis Library 



we mean a collection of 5oth general-purpose mathematical/statistical 
routines and special-purpose data-analysis/pattern recognition 
routines. 

The Image Analysis Library has been partitioned, functionally, 
into the four libraries (and further into packages) as outlined in 
Appendix B. Identification forms were completed for various computer 
subroutines and programs and subsequently arranged by Package (See 
Appendix B) . These completed forms represent an initial attempt 
at defining those routines which will eventually comprise an Image 
Analysis Library for the PDP 11/45.. Care has been taken to insure 
that the identified programs represent what is generally .thought of 
as the best available technique for performing a particular function. 
In a number of cases, several theoretically equivalent computer 
programs yjeve available with minor to extensive variations 
in code. We have exercised our technical judgement in an attempt to 
minimize duplication among the programs appearing in Appendix B. 

In some cases a particular functional requirement appears for 
which no computer program has yet been identified. In those cases 
we have tried to indicate where such programs can be found and 
identified by the individuals who have either developed or are 
developing the programs. 

Each subroutine or program which has been identified is assigned 
an implementation priority according to the following priority 


definitions: 


Priority I - Basic utility subroutine or program 
required for the implementation of 
existing image analysis techniques 
Priority II - Existing image analysis techniques 

(requires majority of Priority I items) 
Priority III - Subroutine or program required to 

support development and/or refinement 
of image analysis techniques 
Priority IV t Image analysis techniques in research and 


development stage. 
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COMPLETED SECTIONS OF PROCEDURES MANUAL 
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PROCEDURES MANUAL OUTLINE 


A. Hardware Documentation 

1. Hardware configuration 
•2. Hardware description by unit 

3. Updating standards 

B. Maintenance. 

1. System log standards 

2. System manager's guide 
3* Hardware 

a. Installation 

b. Check list 

c. Preventive maintenance 

d. Failure 

(1) failure report forms 

(2) check list of repairs in progress 

4. Software 

a. System software 

(1) software performance summary 

(2) software problem report form 

(3) filing and updating standards 

b. Batch and conversational library 

(1) software performance summary 

(2) software problem report form 

(3) filing and updating standards 

c. Interactive software 

(1) . software performance summary 

(2) software problem report form 

(3) filing and updating standards 

d. Software failure bulletin board 
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C. System Software Documentation 

1. Overall system concept 

2 . System programs 

a. DOS Monitor 

b. BATCH - 11 

c. FORTRAN IV 

d. LINK - 11 Linker 

e. LIBR - 11 Librarian 

f. PIP File Utility Package 

g. MACRO - 11 Assembler 

h. EDIT - 11 Text Editor 

i. ODT - HR Debugging Program 

j . Others 

3. Cookbooks 

a. Job control language - BATCH 

b. IMAGE 100 orientation 

4. FORTRAN conversion problems 

5. Concise summary of all error messages issued by system programs 

6. Other system software 

7. Updating standards 

D. Library Documentation 

1. Program documentation standards 

a. Batch and conversational 

b. Interactive special functions 
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2. Image Analysis Library 

a. Pre-analysis data processing library 

b. General mathematical library 

c. General statistical library 

d. Pattern recognition library 

3. Interactive special functions 

4. Updating standards 

E. Interface 

1. Interface documentation standards 

2. Interface documentation by hardware unit 

3. Updating standards 
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COMMENTS ON PROCEDURES MANUAL OUTLINE 


COMMENTS ON PROCEDURES MANUAL OUTLINE 


II 







In the development of this procedures manual outline, we tried to be 
as concise as possible. We make the following comments to guide efforts 
to supply the necessary information. In some instances, we were able to 
identify manuals which partly supply this information. Additional manua ls 
and information, particularly relating to hardware maintenance and system 
software documentation, exist and need to be identified and adapted to fit 
the existing system configuration. 

Documentation of hardware appears in the outline under major divisions 
A, 6 and E, as follows: 

under A (Hardware Documentation), for each unit 

1. General description 

2. Detailed physical description 

3. Power requirements 

4. Operating environment 

5. Normal operations 

6. Principles of operation 

7 • Block diagrams 

8. Electrical diagrams 

9. Accessories 
under B (Maintenance) 

1. Unpacking and incoming inspection 

2. Installation 

3. Warranty information 

4. Initial checkout procedures 
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5. Tools and test equipment required 

6. Preventive maintenance (daily, weekly, quarterly, annually) 


under E 

(Interface) 

1. 

Initialization 

2. 

Interface design 

3. 

Interfacing procedures 

4. 

Timing 

5. 

Data flow 

6 . 

Control panel interface 

7. 

Data lines 

8. 

Code charts 

9. 

Diagnostics 

10. 

Hardware dependent software 


Documentation of programs appears in the outline under major divisions 
B, C and D. By "program” we refer throughout this outline to all types of 
subroutines, functions, programs, program packages, utilities, compilers, 
subsystems and systems. Maintenance of programs is described in B, system 
programs are documented in C (System Software Documentation) and user 
created (library) programs are documented in D (Library Documentation) . 

Software other than programs or program documentation is represented 
throughout the outline. Hardware operations manuals appear in A, diagnostic 
software is documented in E, system operation and management information 
appears in B and initialization procedures for system programs and other 
software procedures are described in C. 

Procedures manual management information is spread throughout the 
manual, mostly in the form of documentation and updating standards. 
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While a reference copy of the entire procedures manual should be kept on 
file , it is clear that some elements need to be located where they can be 
more easily referenced. Documentation standards for the use and maintenance 
of a system library are set in D. 

A. System Hardware Documentation 

1. Hardware configuration (general description) 

a. Overview of hardware with floor plan 

b. Processor (including size of internal storage, cycle time) 

c. Storage media (tape units, disk units) 

d. Output devices 

e. Input devices 

f. Special Image Processing Hardware 

2. Hardware description by unit 

3. Updating standards 

S. Maintenance 

1. System log standards 

These standards specify what is to be kept in the system log, 
the frequency of entries, and those personnel required to 
make and review entries. 

2. System manager’s guide 

DEC manual DEC-ll-OMGRA-A-D (BATCH-11/D0S-11 System Manager’s 
Guide) provides information needed to install and manage 
the operating system. Modifications in the operating system 
necessitated by special interactive requirements must be 
documented here. Non-standard system modules (device drivers 
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and EMT’s) should be described here. 

3. Hardware (by unit) 

a. Installation 

This section includes information on the installation and 
acceptance of the hardware unit; also included is informa- 
tion on warranty, service contracts and general time-limited service. 

b. Checklist 

By this we intend a quick check of the status of the system 
hardware which can be gone over (say) with each change of 
shifts. 

c. Preventive maintenance 

A daily, weekly and so on schedule of preventive maintenance 
should be established and followed by operations personnel. 

d. Failure 

(1) failure report forms 

The failure report procedures currently in use appear 
to be based on those developed for a much larger and 
more general maintenance problem. It appears that time 
could be saved and quality control improved if a modi- 
fication of the current form were made to accommodate 
the special needs of the Image Processing Lab. One 
form, the WITHOLD tag (Form 1055) , seems to be quite 
effective. 

(2) check list of repairs in progress 

This list, now kept in the form of a bulletin board, 
should include parts on order, special equipment 
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or personnel required, the interim status of the 
system and so on. 

4. Software 

The DEC document entitled Software Performance Summary 
contains an extremely attractive and easily maintained 
system for reporting, correcting, filing and updating soft- 
ware problems. We strongly recommend their method replace 
the current method in system software problem reporting, 
and that an adaption to the needs of the Image Analysis 
Library interactive software be made. 

a. System software 

(1) software performance summary 

(2) software problem report form 

(3) filing and updating standards 

b. Batch and conversational library 

(problems in the Image Analysis Library and associated 
documentation) 

(1) software performance summary 

(2) software problem report form 

(3) filing and updating standards 

c. Interactive software 

(problems in the interactive software and associated 
documentation) 

(1) software performance summary 

(2) software problem report form 

(3) filing and updating standards 
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d. Software failure bulletin board 

Here we suggest the establishment of an informal 
bulletin board to alert all system users of software 
failures which have been noted. When a way to avoid 
the problem is know it may be included. Although 
this information may also be included in a.-c. above, 
it seems unlikely most system users would keep up with 
that rather formal documentation. The minor duplication 
of effort involved is negligible compared to the 
benefits. 

C. System software documentation 

As we view the present system configuration, essentially all system 
software is supplied by DEC. We identify here major components of 
system software and make a few suggestions for easier use of the 
system. We also suggest an Image 100 orientation manual. 

1. Overall system concept 

(A general description of system architecture) 

2. System programs 

In each of the system programs mentioned here, we have been 
able to identify a programmer’s guide (a DEC document) but not 
the actual program documentation. For example, the FORTRAN IV 
compiler programmer’s manual contains almost no information on 
how the compiler itself operates. The Procedures Manual 
should contain such material when it is available, 
a. DOS Monitor (Disk Operating System Monitor Programmer's 
Handbook, DEC-ll-OMONA-A-D) 
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b. BATCH-11 (BATCH-11 Users Guide, DEC-ll-OBUDA-B-D) 

c. FORTRAN IV (PDP-11 FORTRAN IV Compiler and Object Time 
System Programmer’s Manual, DEC-ll-LFIVA-A-D) 

d. LINK- 11 Linker, and 

e. LIBR-11 Librarian (BATCH-11/D0S-11 Linker (LINK 11) and 
Librarian (LIBR-11), DEC-ll-ULLMA-A-D) 

f. PIP File Utility Package (BATCH-ll/DOS-11 File Utility 
Package (PIP>, DEC-ll-UPUPA-B-D) 

g. MACRO-11 Assembler (BATCH-ll/DOS-11 Assembler (MACRO-11), 
DEC-ll-OMACA-B-D) 

h. EDIT-11 Text Editor (BATCH-ll/DOS-11 Edit-11 Text Editor 
Programmer’s Manual, DEC- 11-EEDA-D) 

i. ODT-llR Debugging Program (BATCH-ll/DOS-11 ODT-llR 
Debugging Program Programmers Manual, DEC-11-00DA-D) 

j . Others (Other system programs which are operational; it is 
clear this part of the outline is open for expansion as the 
system is added to.) 

Cookbooks 

a. Job control language - BATCH 

By this we intend a folding pamphlet containing basic job 
control statements to compile, link and execute a FORTRAN IV 
program in the BATCH mode. As actual operating experience 
is gained with the BATCH system, it seems likely that revision 
of the pamphlet will be necessary. 

b. IMAGE 100 orientation 

A brief applications oriented description of use of the 



IMAGE 100 system functions in actual application with color 
photographs of how the screen looks, examples of what choice 
of channels or passes is appropriate for various users; this 
information would he of great value to the operators in ex- 
plaining to a new user what to expect, and would also serve 
to illustrate to potential users some of the actual applications. 

4. FORTRAN conversion problems 

The conversion of a program from one computer installation to 
another is never trivial, but, when the conversion involves 
change to a 16 bit machine from a machine with at least 32 bit words, 
the problems may seem almost insurmountable. This section of the 
outline should initially contain hints on probable problems in 
converting an operational FORTRAN IV program to the PDP 11/45 
computer. In order to make this information as useful as possible, 
it should be made easy to revise and add to as experience is 
gained with actual conversion problems. 

5. Concise summary of all error messages issued by system programs 
A tabular summary of error messages enabling a user to discover 
which component of the system issued the message and its meaning 
should be available to both BATCH and Interactive users. 

6. Other system software 

7 . Updating standards 

D. Library Documentation 

Since all programs must be tested for correctness and occasionally 
modified, adequate documentation is a must so that future users may know 
exactly what a particular program is supposed to accomplish and how. 
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1. Program documentation standjirds 

a. Batch and conversational 

b. Interactive special functions 

2. Image Analysis Library 

a. Pre-analysis data processing library 

b. General mathematical library 

c. General statistical library 

d. Pattern recognition library 

3. Interactive special functions 

These are interactive programs run from the IMAGE 100 console. 

4. Updating standards 

E. Interface 

In a PDP 11/45 computer system, hardware units generally communicate 
with the Central Processing Unit and each other through the Unibus high speed 
data bus by direct interfacing (wires) or by data communication interfacing 
(usually a telephone or other low-speed line with a moden unit on each end) . 
Each unit of hardware communicates with an interactive user or with the 
computer operator through a control panel interface. Especially in an 
interactive environment, programmers (and sometimes users) need to be 
aware of technical details of how commands are decoded, how data is routed 
and how correct timing is maintained. Subroutines and functions exist which 
would allow a FORTRAN IV programmer to write interactive programs making 
efficient use of interactive hardware capabilities if the progr amm er' 
understood how to use the interfacing software. 

1. Interface documentation standards 

Each type of interface should be provided with a documentation 
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standard sufficient to allow the interface to be maintained 
and used. 

2. Interface documentation by hardware unit 

3. Updating standards 
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COMPLETED SECTIONS OF PROCEDURES MANUAL 


A. 3 


B.4.a (3) 
B.4.b (2) 
B.4.b (3) 
B.4.c (2) 
B.4.C (3) 


C.3.a 


C.4 


C.7 


D.l 


D.4 


E.l 


E.3 


Hardware Documentation Updating Standards 

System Software Maintenance — Filing and Updating 
Standards 

Batch and Conversational Library Maintenance — 
Software Problem Report Form 

Batch and Conversational Library Maintenance — 

Filing and Updating Standards 

Interactive Software Maintenance — Problem Report 
Form 

Interactive Software Maintenance — Filing and Updating 
Standards 

Job Control Language Cookbook 
FORTRAN Conversion Problems 

System Software Documentation Updating Standards 
Library Documentation Standards 
Library Documentation Updating Standards 
Interface Documentation Standards 
Interface Updating Standards 
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A. 3 Hardware Documentation Updating Standards 

Each change in the system hardware configuration shall be documented 
as follows: 


Change 

Action 

Addition of a 
hardware unit 

Appropriate information shall be added to the 
Procedures Manual. (For reference, the distri- 
bution of this information throughout is listed 
below.) Section A.l (Hardware configuration) 
shall contain references to sections of the 
Procedures Manual containing information about 
the added unit immediately following the 
general description of the unit. 

Removal of a 
hardware unit 

! 

No longer applicable information shall be 
indicated by means of vertical line along the 
outer margin with a dated comment describing 
the change as "removed (date)". When several 
consecutive pages which would require this 
treatment can be fastened together securely 
only the first page need have the change 
indicated. 

Modification or 
rearrangement 

Changes induced shall be indicated by means 
of a vertical line along the outer margin 
with a dated comment describing the change as 
"changed (date) (page reference)" with a 
reference to currently applicable information. 
When possible, additions shall be physically 
near (in the Procedures Manual) the items they 
modify. 
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Reference information on distribution of hardware documentation: 


Item ' Section 


Accessories 

A. 2 

Block diagrams 

A. 2 

Code charts 

E.2 

Control operation 

A. 2 

Control panel interface 

E.2 

Data flow 

E.2 

Data lines 

E.2 

Detailed physical description 

A. 2 

Diagnostic software 

E.2 

Electrical diagrams 

A. 2 

Floor plan 

A.l 

General description 

A.l 

Hardware dependent software 

E.2 

Initial checkout procedures 

B . 2. a 

Initialization (interface) 

E.2 

Installation 

B.2.a 

Interface design 

E.2 

Interfacing procedures 

E.2 

Logic flow 

A. 2 

Normal operation 

A. 2 

Operating environment (physical) 

A.l 

Power requirements 

A.l 

Preventive maintenance 

B,3.c 

Principles of operation 

A. 2 

Timing 

E.2 

Tools and test equipment required 

B.3.c 

Unpacking and inspection 

B.3.a 

"Warranty information 

B.3.a 
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B. 4. a (3) System Software Maintenance-Filing and Updating Standards 

Filing and updating standards for system software failure report 
forms shall he those described in DEC document Software Performance 
Summary (DEC-ll-XSPSB-A-D) , Copyright 1974, Digital Equipment Corporation, 
§ 2 . 
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B.4.b (2) Batch and Conversational Library Ms in t enanc e — Software Problem 
Report Form 

The maintenance of the Image Analysis Library requires the design 
and use of a new program problem report form, patterned after the 
system used by DEC (appearing in the Procedures Manual in B.4.a (2)). 

Each problem report takes the following overall form; each part will be 
explained in detail. 


— 

DATE 

I. TITLE 


II. SUBTITLE 


III. PROBLEM: 

USER: 

IV. SOLUTION: 

PROGRAMMER: 


V. CODING 
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DATE: The date of the current version, of the problem report. 

I. TITLE 

The title of the problem report refers to the component of the 
Image Analysis Library involved (for example, the Pattern 
Recognition Library) 

XI. SUBTITLE 

This brief statement (which should fit on one line) , gives the 
user a hint about the content of the problem report. The 
title and subtitle are to be used in a table of contents. 

III. PROBLEM USER: 

A paragraph is used to describe the problem in general terms, 
perhaps with an example, warnings, etc. The name and organiza- 
tion of the user reporting the problem appear here. 

IV, SOLUTION: PROGRAMMER: 

If there is a way to avoid the problem it is stated here. 

V. CODING 

This block is used for filing purposes. Below is a close-up 
view of the coding block 


PROGRAM NAME 

PROGRAM NUMBER 

SEQUENCE # 

ADDITIONAL INFORMATION (SUBROUTINES, FUNCTION) 

PAGE 

OF 

NEW REPLACEMENT H 

ORIGINAL DATE 



......... 
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Problem reports are filed by program (subroutine, function) name. 
Each program in the Image Analysis Library is assigned a number when 
accepted. The sequence numbers are assigned by program problem; as 
each new problem report for that program is filed it is assigned the 
next highest sequence number . Any additional information which may 
help a user is recorded. If the problem report is new (for this problem) 
an ’x' is placed in the NEW block; if a replacement, the revision 
number is placed in the REPLACEMENT # box. The original date is the 
date of the first report on this problem. 
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B.4.b (3) Batch and Conversational Library Maintenance — Filing and 
Updating Standards 

Program Problem Reports shall be filed by program name and indexed 
by function as an element in the Image Analysis Library. When a problem 
is fixed or the program ■withdrawn, or information is added, the problem 
report shall be removed or replaced by a revision. 
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B.4.c (2) Interactive Software Maintenance — Problem Report Form 

Interactive use of the system is by users who may not be programmers. 
To account for this, a slight redesign of the form in B.4.b (2) has been 
made: 


I. 

MENU 

DATE 

II. 

OPTION/FUNCTION 


III. 

SUBTITLE 


IV. 

PROBLEM 

USER: 

i V. 

SOLUTION 

PROGRAMMER: 



VI. CODING 
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DATE: Date of the current version of the problem report. 

X. MENU 

The menu refers to a user-selected set of options grouped 
together by the interactive system. 

XI. OPTION/ FUNCTION 

This describes which option was selected of the list under the 
menu of (I) above. 

III. SUBTITLE 

This is a brief (one line) description of the problem. 

IV. PROBLEM: USER: 

This is a one paragraph description of the problem with an 
example of a failure followed by the name and organization; action 
of the user reporting the problem. 

V. SOLUTION: PROGRAMMER: 

If there is a way to avoid the problem it is stated here, followed 
bv the name and organization of the programmer who devised the 
solution. 

VI. CODING 

This block is used for filing purposes. Below is a close-up view of 
the coding block. 


OPTION 

SEQUENCE # 

ADDITIONAL INFORMATION 

PAGE 

OF 

NEW REPLACEMENT # 

ORIGINAL DATE 
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Problem reports are filed by option name. Sequence numbers are 
assigned by option; as a new problem report for an option is reported, 
the next highest sequence number is assigned. If the problem report is 
new, an *x* is placed in the NEW block; otherwise the revision number 
is placed in the REPLACEMENT # box. 
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B.4.c (3) Interactive Software Maintenance — Filing and Updating Standard 


Problem reports shall be filed by option name and indexed by menu 
and function as an Interactive special function. When a problem is fixed 
or more information is added, the problem report Bhall be removed or 
replaced by a revision. 
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C.3.a Job Control Language Cookbook 
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I 

C.4 FORTRAN Coaversion Problems j 

I 

The principal difficulties which arise when conversion of a FORTRAN 

; 

IV program from a different implementation to PDP 11/45 is attempted 

stem from word length differences, core size differences, JCL dif- ■ 

ferences and actual compiler bugs. The purpose of this section of 
the Procedures Manual is to allow users to share experience gained in 

i" 

attempting conversion so that an overall saving in effort (and j 

frustration) can be realized. Much of the following information is 
contained In the FORTRAN Manual and the Software Performance Su mma ry . 

INTEGER Size: 

Since the word length of the PDP 11/45 is 16 bits, 
the maximum signed integer is 2**15 - 1 = 32,767. 

"Double Precision" integers are not supported. 

BYTE (LOGICAL * 1) Variables and Arrays: 

BYTE (or LOGICAL * 1) variables and arrays cannot be 
initialized in DATA statements. The compiler allows 
simple BYTE variables to be used in arithmetic 
statements but does not allow them in arithmetic IF 
statements. Comparison with Hollerith strings leads 

to unexpected results because Hollerith strings are j 

treated as type INTEGER. A BYTE array element used as \ 

a subroutine argument sometimes causes an execution 
time addressing error. 

ENCODE /DECODE Problems: 

Core-to-core transfer under FORMAT control require care 
in that the area of core must be contiguous. Thus if 
INTEGER arrays are being processed, the program will have 
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to have been compiled with the /ON switch. 

F FORMAT Problem: 

Double precision variables are truncated when output 
under F FORMAT control. 

EQUIVALENCE Problem: 

Alignment to even bytes is necessary in EQUIVALENCE statements. 
DEFINE FILE Problem: 

Multiple DEFINE FILE statements are not compiled correctly, 
nor is a diagnostic message issued. 

UNARY MINUS: 

Unary minus is given precedence over exponentiation. Thus 
-2**2 is treated as (-2)**2 instead of the standard -(2**2). 
ENDFILE Problems: 

ENDFILE cannot be the statement to be conditionally 
executed in a logical IF. ENDFILE cannot be the first 
executable statement in a subroutine. 

Logical Unit Numbers as Subroutine Arguments: 

The compiler does not- Recognize integer variables passed 
as subroutine arguments as valid I/O logical unit numbers. 

A dummy assignment to a local variable is necessary to fix the 
problem. 
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C.7 System Software Documentation Updating Standards 


As system software is updated or revised, changes shall be indicated 
in the Procedures Manual as follows: 


I 


Change 


Action 


Addition of a 
component of 
system software 


Removal of a 
component of 
system soft- 
ware 


Modification 


Documentation and user information shall be 
included at the proper place. 


No longer applicable documentation shall be 
indicated by a vertical line along the outer 
margin with a dated comment describing the 
change as "removed (date)". When several con- 
secutive pages which would require this 
treatment can be fastened together securely 
only the first page need have the change 
indicated . 


Changes induced shall be indicated by a 
vertical line along the outer margin with 
appropriate dated comment and reference to 
corrected or updated documentation. 



D.l Library Documentation Standards 

The sharing of computer programs by a number of users requires that 
their documentation be managed so as to insure uniformity without imposing 
an undue clerical burden. While documentation must be performed at the 
same time the program is being developed, a requirement that a formal 
standard be satisfied at all levels would be costly. We recognise 
four levels of documentation: 

LEVEL 1 Minimal level 

Level 1 documentation standards apply to single-use programs of minimal 
complexity. The documentation required shall be the completion of the 
first page of the computer program documentation form, a liberally 
commented listing and test case. 

LEVEL 2 Internal level 

Level 2 documentation standards apply to larger special purpose programs 
which appear to have no sharing potential. The documentation required 
shall be the completion of pages 1 - 4 of the computer program documentation 
form, a liberally commented listing and a test case. 

LEVEL 3 Working document level 

Level 3 documentation standards apply to programs which are intended 
to be shared, either as parts of the Image Analysis Library or as 
special functions for Interactive use. In addition to Level 2 
documentation, all basic elements of good documentation must be 
prepared and kept in the Library. 
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LEVEL 4 Interactive use level 


Level 4 documentation standards apply to programs which are intended 
to be used in an interactive mode. In addition to Level 3 documentation, 
page 5 of the computer program documentation form (which describes 
operation of the program in an interactive mode to a user who may not under- 
stand the program) must be completed. 

Basic elements of good documentation: 

Title page 
Table of contents 

Abstract - brief summary of function or purpose 

Problem task description - technical description of problem to be solved 
or task to be accomplished 

Method of solution - mathematical techniques used, logical sequencing 
within the program and data files created by the program 
Program description: 

Operating Environment - hardware and software requirements 
Program specifications - a detailed description of the programming 
techniques used in writing the program; e.g. calling sequence, 
overlay structure, test plan, COMMON usage 
Subprograms (other than library) - with appropriate level of 
do cumentatxon 

Non ANSI FORTRAN IV library subprograms — names and brief 
description 
Source code listing 
Flow charts — general and detailed 
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Operating instructions - description of procedure required to 
set up and run the program, including: 

Deck setup - showing placement of control cards, source or 
object decks and data 

Description of input - file, record and data element 

descriptions and formats including the origin of each 
data element 

Description of output - file, record and data element 

description and formats of each data element including 
scratch or intermediate files 

Restrictions or limitations - hardware or software restrictions, 
data ranges and capacities, program behavior when 
restrictions are violated 

Accuracy characteristics - including dependence on machine 
word length or floating point arithmetic 

Diagnostics emitted - including probable cause of program- 
generated messages 

Test Cases 

Special instructions - interactive mode 

Symbols - a complete tabular list of all symbols and array 
names used in the program, defined with reference to 
the mathematical or technical writing notations and terms 
used in the problem description, with units given where 
applicable 

Record of changes - arranged in columnar form providing spaces for 
the change number, date and individual making the change. 


PM-33 


with a brief description of the reason for the change; 
actual changes should be made in the document itself 
References 
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Computer Program Documentation Form 


Date: / / Program No. 

Parent Program No. 

Keywords and phrases: 

Title: 

Author: 

Author location: 


Contact : 



Application: _ 

Computer type: System or monitor: 

Source language(s) : _( %) ( %) 

Memory required, other than monitor: 

List program numbers of any required subprograms that are documented separately 
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References 


Revision or modification history: 


Date 


Reason 
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" " "I 

Interactive user instructions j 

Menu j 

l 

Option name 

Load and start procedure 
Expected wait time 
User instructions 


Error messages l 

I 

Error recovery j 

i 

i 


Effect on system parameters 

! 

UIC j 

! Protection code 

i 

i 

I 

Files created j 

COMMON variables . j 

Image 100 hardware affected 
Other 
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D.4 Library Documentation Updating Standards 

Each program in the Image Analysis Library and each Interactive 
special function shall be documented with a revision or modification 
history in which a record is kept of the revision number, date, 
programmer and reason. Should a program which is on a public 
disk be withdrawn, notice shall be posted on the software failure 
bulletin board. Every possible step shall be taken to insure 
upward compatibility of program revisions or modifications. 
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E.l Interface Documentation Standards 


a. Hardware - CPU interface 

The connection of hardware units to the central processing 
unit shall be described in detail including information on 
initialization, data flow, data lines, timing, interface design 
and interfacing procedures. 

b. Hardware - control panel interface 

The connection of control panel switches, keys and controls to 
the data flow and interface bit status shall be described in 
detail. 

c. Hardware - software interface 

In the documentation of each Interactive software unit, the 
action of hardware controls on software variables shall be 
described. Utility software which exists for the purpose of 
enabling users to poll the status of hardware units or hardware 
controls shall be documented fully here. 


G 


E.3 Interface Updating Standards 

The updating standard for change in the hardware configuration is set 
in section A. 3. 

Changes in the utility software shall be documented following the 
standard for system software maintenance under the component 
name "UTILITY". 
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APPENDIX B 

IMAGE ANALYSIS LIBRARY OUTLINE 
COMPUTER PROGRAM IDENTIFICATION FORMS 
IMPLEMENTATION PRIORITY SUMMARY 
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IMAGE ANALYSIS LIBRARY OUTLINE 


PRE-ANALYSIS DATA PROCESSING LIBRARY 

1. PREPROCESSING PACKAGE 

a. RADIOMETRIC CORRECTION 

b. SUN ANGLE CORRECTION 

c. HAZE CORRECTION 

d. MEAN LEVEL ADJUSTMENT 

e. GENERAL SCANNER UNIT CORRECTION (IMAGE 100) 

f. ERTS DISTORTION CORRECTION 

g. I/O OF VARIOUS FORMATS 

2. REGISTRATION PACKAGE 

a. SINGLE PASS TO GEOGRAPHICAL MAP OR AERIAL PHOTO 

b. SEVERAL PASSES TO DESIGNATED PASS 

3. TRAINING FIELD SELECTION PACKAGE 

a. IRREGULARLY SHAPED FIELDS 

b. HISTOGRAMS FROM LINEAR COMBINATIONS OF n CHANNELS 

c. GREY SCALE MAP PROGRAM 

d. CONTOUR PROGRAM 

4. SIMULATED DATA GENERATOR PACKAGE 

a. RANDOM NUMBER GENERATOR-VARIOUS DISTRIBUTIONS 


b. SIMULATED ERTS DATA 
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IMAGE ANALYSIS LIBRARY OUTLINE (CONT'D) 

B. GENERAL MATHEMATICAL LIBRARY 

1. LINEAR ALGEBRA PACKAGE 

a. LU DECOMPOSITION 

b. MATRIX INVERSE 

c. CHOLESKY METHOD' — SYMMETRIC POSITIVE DEFINITE 

d. GENERALIZED INVERSE 

e. EIGENVALUE PROBLEM: Ax = Ax, A SYMMETRIC 

f. GENERALIZED EIGENVALUE PROBLEM: Ax = ABx, A, B SYMMETRIC 

g. HOUSEHOLDER TRANSFORMATION 

h. SINGULAR VALUED DECOMPOSITION 

i. POLAR DECOMPOSITION 

j . MATRIX MANIPULATION ROUTINES (GENERAL) 

2. OPTIMIZATION PACKAGE 

a. LINEAR PROGRAMMING 

b. QUADRATIC PROGRAMMING 

c. GRADIENT-STEEPEST DESCENT 

d. FLETCHER POWELL DEFLECTED GRADIENT 

e. NEWTON’S METHOD AND VARIATIONS 

f. CONSTRAINED LEAST SQUARES 

3 . APPROXIMATION PACKAGE 

a. SPLINE APPROXIMATION — S INGLE VARIABLE 

b. SPLINE APPROXIMATION— TWO VARIABLES 

4. TRANSFORM PACKAGE 

a. FAST FOURIER TRANSFORM 

b. HADAMARD TRANSFORM 

c. WALSH TRANSFORM 
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IMAGE ANALYSIS LIBRARY OUTLINE (CONT'D) 


C. GENERAL STATISTICAL LIBRARY 

1. STATISTICAL SUMMARY PACKAGE 

a. MEAN VECTOR COMPUTATION 

b. VARIANCE-COVARIANCE MATRIX COMPUTATION 

c. MODE, MEDIANS, RANGE, ETC. 

2. DENSITIES & DISTRIBUTIONS PACKAGE 

a. MULTIVARIATE NORMAL 

b. UNIVARIATE NORMAL 

c. ERROR FUNCTION AND INVERSE ERROR FUNCTION 

d. UNIVARIATE DENSITIES— MISCELLANEOUS 

3. REGRESSION PACKAGE 

a. MULTIPLE LINEAR REGRESSION 

b. NONLINEAR REGRESSION 

4. STATISTICAL TESTS PACKAGE 

a. HYPOTHESIS TESTS 

b. GOODNESS OF FIT 

c. ANALYSIS OF VARIANCE 



IMAGE ANALYSIS LIBRARY OUTLINE (CONT'D) 

D. PATTERN RECOGNITION LIBRARY 

1. CLUSTERING PACKAGE 

a. ISOCLS 

b. SEMI-SUPERVISED CLUSTERING (G.E.) 

2. CLASSIFICATION PACKAGE 

a. BAYES OPTIMAL— m GAUSSIAN CLASSES 

b. BAYES OPTIMAL— CONVEX COMBINATIONS OF GAUSSIAN 

c . TABLE LOOKUP 

d. MAXIMUM LIKELIHOOD RESOLUTION (G.E.) 

3. FEATURE SELECTION PACKAGE 

a. DIVERGENCE CRITERION 

b. BHATTACHARYYA DISTANCE 

c. PROBABILITY OF MISCLASSIFICATION 

d. LINEAR DISCRIMINANT FUNCTIONS 

e. TRANSFORMED DIVERGENCE 

4. ESTIMATION OF PROPORTION PACKAGE 

a. MAXIMUM LIKELIHOOD 

b. METHOD OF MOMENTS 

c. ESTIMATED CONFUSION MATRIX 

d. CLASSIFICATION 

e. BAYES WEIGHTS— CONFUSION MATRIX 

5. SPATIAL PATTERN RECOGNITION PACKAGE 

a. FIELD BOUNDARY DETECTION 

b. SPATIAL FEATURE DEFINITION 

c. CLASSIFICATION USING SPATIAL FEATURES 

6. SIGNATURE EXTENSION PACKAGE 

a. MAXIMUM LIKELIHOOD FOR MEAN, COVARIANCE, PROPORTIONS 
USING UNLABELED SAMPLES AND CONTRACTABILITY 

b. Ax + b, A DIAGONAL, CORRECTION TRANSFORMATION 
(UNLABELED SAMPLES) 
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COMPUTER PROGRAM IDENTIFICATION FORMS 

The computer subroutines and programs which were identified 
appear in the following pages. They have been arranged into the 
packages defined in the outline of Image Analysis Library. In some 
cases no computer programs were Identified which provide a required 
functional capability. We have tried to indicate through notes on 
the header page for each package where such programs exist for purposes 
of future identification and acquisition. 
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A„ PRE-ANALYSIS DATA PROCESSING LIBRARY 
1. PREPROCESSING PACKAGE 


Computer programs are available to meet some of the functional 
requirements of this package and can be supplied by either LEC or ERIM, 
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COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 2 


I. NAME: GENERAL RADIOMETRIC CORRECTION 

II. SOURCE (Include coder/person with knowledge of program usage) : 

General Electric/Jim Brienly 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: IMAGE 100/PDP-ll 

IV. PURPOSE: See next page 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See next page 

Description: See next page 

VI. SPECIAL FEATURES: None 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: X A vailabl e N ot Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: $7,712 

I s x I s not compatible with Image 100 configuration In present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e X Special purpos e X P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: X LACIE X Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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• *1.1.2 GENERAL RADIOMETRIC CORRECTION 

Adjustment of the radiometric amplitudes of raw ERTS-I MSS data may be effected 
using this program. Three modes of operation are available for handling a variety of 
data adjustment techniques: 

a. General Mode 

This mode assumes that the data requii-es correction as a function of x- y 
position and channel. The form of the correction applied to each pixel 
is. 

S' (x,y) = S(x,y)*f(ch)*s(ch)-ft(x,y)*g(ch) 

where 

S is the radiometric amplitude of each pixel; x,y are ERTS MSS frame 
coordinates; ch is the ERTS MSS channel; f, g, s, t are user-supplied 
functions . 

b» Channel Correction Mode 

This mode permits the user to apply a stored lookup table correction 
to each of the four channels. This is to correct ( or linearise) non- 
linear scanner sensor outputs. The lookup tables are generated under 
* user control for each channel from the following functions: 

a) S’ = aS® + bS + c (square law function) 

(2) S’ = a In (S+l) + b (log function) 

(3) S' = ae s + b (exponential function) 

(4) Three-segment piecewise linear curve with break-points at a and b 
corresponding to the 10 percent and 90 percent amplitude points. 

-c. Sensor Correction Mode 

This mode is applied to correct for calibration differences between 
ERTS-1 MSS sensors. Each image line is referenced back to its 
sensor and a stored lookup table correction is applied. The lookup 
tables may be loaded at the user's request; this is accomplished by 
scanning the image and loading each table with the corrections that pro- 
vides uniform values lor the mean of each sensor within a band. In this 
fashion, so-called "banding" corrections may be made. 
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l. PRE-ANALYSIS DATA PROCESSING LIBRARY 
2. REGISTRATION PACKAGE 


A. 2-1 


COMPUTER PROGRAM IDENTIFICATION FORM 
THIS FORM PREPARED BY LEG 


1 of 2 


I. NAME: CORLAT 

II. SOURCE (include coder/person with knowledge of program usage): 

S. YAO/LEC 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: UNIVAC 1108 

IV. PURPOSE: Does two-dimension image correlation using the fast Fourier 
transform technique 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 

Description: By inputting a reference and an overlay image and the 

correlation patch-size, the output correlation function and 
correlations coefficients are obtained. 

VI. SPECIAL FEATURES: One can either correlate on designated Line-pixel 

points or specify a set of correlation grid points. 

VII. SUBROUTINES REQUIRED: HARM (2D-FFT from IBM Scientific Subroutine package) 

VIII. DOCUMENTATION: A vailabl e x N ot Available Will be Available 

IX. _x Software available off-the-shelf. Approximate Cost: 

I s -s- I s not compatible with Image 100 configuration In present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

x O ther (specify) Pre-processing — Registration 
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COMPUTER PROGRAM IDENTIFICATION FORM 


2 of 2 


I. NAME: REGSTR 

II. SOURCE (include coder /person with knowledge of program usage): 

S. YAO 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: UNIVAC 1108 

IV. PURPOSE: Does linear (1st order) geometrical corrections on 2-D images 

based on six input parameters — two scaling, two translation and two 
rotation. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference : 

Description: REGSTR mapd the overlay image onto the reference image, 

with output image similar to the reference image in 
geometrical properties. A buffer is created to store the 
overlay image for output one line at a time. 

VI. SPECIAL FEATURES: REGSTR does "REGISTRATION" based on results obtained 

from CORLAT. It contains the tape-merge provision, 
that is, one output image can be obtained from merging 
two adjacent ERTS strips. 


VII. SUBROUTINES REQUIRED: 


VIII. DOCUMENTATION: Available v Not Available Will be Available 

IX. Software available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

v PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 


XI. 


APPLICATION : LACIE x Non-LACIE 

Training Classification Estimation of Proportion 


Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 


x O ther (specify) pre-processing — REGISTRATION 
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A. PRE-ANALYSIS DATA PROCESSING LIBRARY 


3. TRAINING FIELD SELECTION PACKAGE 


A program which provides the capability of selecting irregularly 
shaped training fields has been developed on the Image 100 System by 
Ted Kell, LEC. 

A display package which makes use of linear combinations of 
channels is being developed by Dr. Jack Bryant, Texas A&M University. 
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COMPUTER PROGRAM IDENTIFICATION FORM 
THIS FORM PREPARED BY LEC 


1 of 1 


I. NAME: IRREGULARLY SHAPED TRAINING FIELDS 

II. SOURCE (Include coder/person with knowledge of program usage) : 
Ted Kell 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: IMAGE 100 

IV. PURPOSE: Select training fields with odd shapes and/or sizes. 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: Program under development by Ted Kell. 


Description: This program will permit the selection of odd shaped training 
fields (the current system permits only rectangular fields). 


VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: Availabl e X N ot Available Will be Available 

IX. Software available off-the-shelf. Approximate Cost: 

I s x I s no t c ompatible with Image 100 configuration in present form 
Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : x LAC IE X Non-LACIE 

x T raining Classification Estimation of Proportion 

F eature Selection Signature Extension 

Temporal Sampling Systems Developmen t X T echnique Development 

Crop Assessment Yield Estimation 

Other ' (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
THIS FORM PROVIDED BY LEC 


1 of 1 


I. NAME: Line find 

II. SOURCE (include coder /person with knowledge of program usage) : 

Ted Kell 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: PDP11/45 

IV. PURPOSE: To detect lines in a scene and display them on the CRT. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 

Description: 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: Available Not Available Will be Available 

IX. Software available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpos e x P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : x LAC IE x Non-IACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 4 


I. NAME: Grey Scale Map Package For High Speed Printer 

II. SOURCE (include coder/person with knowledge of program usage): 

Jack Bryant, Dept. Math. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV, PLI COMPUTER: IBM 360/65 

IV. PURPOSE: To produce single channel grey scale maps on a high speed 
printer. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: Jack Bryant, Program documentation OCM, HISTOG, PRINTUM, 

Report # 35, Contract NAS-9-12777, Department of Mathematics, University 
of Houston, August, 1974. 

Description: Computes histograms from single channel of ERTS data 
(single pass or registered) and produces a grey scale map on the high 
speed printer (8 levels of grey) , 

VI, SPECIAL FEATURES: Three separate main programs - not a subroutine 
package. 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: -0- 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

x I n language other than FORTRAN. Specify IBM 360 Assembly Lang uage 
Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpos e x P rogram package 

PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : x LACIE x Non-LACIE 

x T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 


2 of 4 


I. NAME: PROGRAM PRINTUM 

II. SOURCE (Include coder /person with knowledge of program usage) : 

Jack Bryant, Dept. Math. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360/65 

IV. PURPOSE: Calculate cumulative distributions of the histograms. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: Jack Bryant, "Program Documentation OCM, HISTOG, PRINTUM," 
Univ. of Houston Report 35, August 1974, NAS9-12777. 

Description: Calculate cumulative distributions (estimate maximum & 
minimum significant data values to be input to OCM) 

VI. SPECIAL FEATURES: , 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION : x A vailable Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: -0- 

I s xl s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

x T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Samplin g S ystems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 


3 of 4 


I. NAME: PROGRAM HISTOG 

II. SOURCE (include coder/person with knowledge of program usage) : 

Jack Bryant, Dept. Math. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 3^/65 

IV. PURPOSE: Outputs complete histograms of all channels. 

V. ALGORITHM/TECHNIQUE (short description; Include reference if available) 

Reference: Jack Bryant, "Program Documentation OCM, HISTOG, PRINTUM’.’ , 
Univ. of Houston Report 35, Aug. 1974, NAS9-12777. 

Description: Accumulates, prints jpunches, complete histograms of all 
channels used as input to OCM. 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General ( rpos e x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

jc Training Classification Estimation of Proportion 

Feature Selection Signature El tension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: PROGRAM OCM 

II. SOURCE (Include coder/peraon with knowledge of program usage): 

Jack Bryant, Dept of Math. , TX A&M Unlv 77843 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360/65 

IV. PURPOSE: One channel map. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: Jack Bryant, "Program Documentation OCM, HISTOG, PRINTUM", 
Unlv. of Houston Report 35, Aug. 1974, NAS9-12777. 

Description: Reproduces digital picture on high speed printer by printing 
an appropriate selection of characters. 

VI. SPECIAL FEATURES: 8 shades of gray can be reproduced. 

Requires high speed printer 

VII. SUBROUTINES REQUIRED :HP05, COLS, TODEC 

. VIII. DOCUMENTATION: x A vailable Not Availabl e W ill be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

x S ome machine language conversion necessary 

x In language other than FORTRAN. Specify IBM 3f0 Assembly Langu age 
Software must be designed and developed for Image 100 

IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

APPLICATION: x LAC IE x Non-LACIE 

x T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) _ 


X. 


XI. 
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I. 

II. 

III. 

IV. 

V. 


VI. 

VII. 

VIII. 

IX. 


X. 

XI. 


NAME: K.ONTUR Program Package 

SOURCE (include coder/person with knowledge of program usage) : 

Charles Eleuterius, Dept, of Oceanography . TAMU 
PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

PURPOSE: To plot, on either a CALCOMP or Gerber plotter, a 
contour map built from a two-dimensional array of data 
values. The array size may be quite large. 

ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: KONTUR, a modification to program CONMAP by 
David C. Haley of Lockheed Electronics Company written 

under NASA contract #NSA 0-5384. For reference to modification — "source" above. 
Ueswription: The surface is approximated by elementary rectangles with 
grid points at the corners. All isolines that traverse the rectangle 
will be drawn while those four points are being considered. 


SPECIAL FEATURES: A code deck is read as the first part of the input 
data; it contains a decision table to determine cases. The code 
deck must be obtained with the program. 


SUBROUTINES REQUIRED: CONMAP, READ IN, TAMU 


DOCUMENTATION: x A vailable Not Availabl e W ill be Available 

x S oftware available off-the-shelf. Approximate Cost: -7- 

Is Is not x c ompatible with Image 100 configuration In present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

IMPLEMENTATION LEVEL 

x G eneral purpose Special purpos e x P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

APPLICATION: x LAC.IE y Non-LACIE 

x T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development x Technique Development 

Crop Assessment Yield Estimation 

Other (specify) . __ 
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I. NAME: MAIN 

II. SOURCE (Include coder /person with knowledge of program usage): 

Charles Eleuterius 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Inputs code deck and control parameters, calls CONMAP. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See KuNTUR reference. 

Description: 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: CONMAP, LINE4 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: _________ 

Is Is no t x compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : X LAC I E X N on- LAC IE 

x T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE CONMAP 

II. SOURCE (include coder /person with knowledge of program usage) : 

Charles Eleuterius 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To permit contouring using very large data arrays. The output 
contour maps will be drawn on the CALCOMP or Gerber plotter. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See KONTUR reference 

Description: Four data points are considered at a time, connected to 
form two triangles. The intersection line of these triangular planes 
with each of the contour planes is plotted. 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED : READ IN, PLOT, FACTOR, SYMBOL, TAMU 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is not x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non- LAC IE 

x T ra in ing Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE READ IN 

II. SOURCE (Include coder/person with knowledge of program usage): 

Charles Eleuterius 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 1160/65 

IV. PURPOSE: To input grid array data 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See KONTUR reference 

Description; 

VI. SPECIAL FEATURES: User written so that data may be entered in form to 
conform to and take advantage of plotter characteristics of travel 
direction and resolution. 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

x T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t y T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE TAMU 

II. SOURCE (include coder /per son with knowledge of program usage) : 

Charles Eleuterius 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To consider each successive rectangle and determine the type 
of isoline case. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See KONTUR reference. 

Description: 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: DRAW 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost; 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

x T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE DRAW 

II. SOURCE (Include coder/person with knowledge of program usage) : 

Charles Eleuterius 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To draw actual contours, labeling the isolines, for each 
rectangle. 

V. ALGORITHM/ TECHNIQUE (short description; Include reference if available) 
Reference: See KONTUR reference. 

Description: 

VI. SPECIAL FEATURES: 

* 

VII. SUBROUTINES REQUIRED: PLOT 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX, x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: xLACIE x Non- LAC IE 

x T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Samplin g S ystems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE MINVSP (A, N, E, K) (single precision) 

SUBROUTINE MINVDP (A, N, E, K.) (double precision) 

II. SOURCE (include coder /person with knowledge of program usage) : 

L. F. Guseman, Jr. , Math. Dept. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360/65 

IV. PURPOSE: To invert a nonsingular N*N matrix A 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: E. G. N. Clayton, Compact methods for inverting matrices and 
solving simultaneous equations by use of Gauss-Jordan elimination, 

M. S. Thesis, Dept. Math., TAMU, May, 1962. 

Description: 

Uses Gaussian elimination with pivoting. 

VI. SPECIAL FEATURES: 

The inverse is built into the array originally occupied by 
A. The original matrix A is lost. 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: -?- 

I s x Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x ppp 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

T rainin g x C lassificatio n x E stimation of Proportion 
x F eature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t y T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: LU-Decomposition Package 

II. SOURCE (include coder/person with knowledge of program usage): 

N.W. Naugle, Math. Dept. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360/65 

IV. PURPOSE: Package for solving linear systems AX=B where A is a square 
nonsingular matrix. Also provides capability to compute A“^ 

V, ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: George Forsythe and Cleve B. Moler, Computer Solution 
of Linear Algebraic Systems, Prentice- Hall, Inc. Englewood 

i 1967 > w- 68 " 72 - 

Computes LU-decomposition of A 

VI. SPECIAL FEATURES: 

Original A is saved. Iterative improvement of solution or inverse 
is built in. 

VII. SUBROUTINES REQUIRED: 

DECOMP, SOLVE, IMPRUV, SING 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware, available off-the-shelf. Approximate Cost: -?- 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpos e x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : x LACI E x N on- LAC IE 

Trainin g x C lassificatio n x E stimation of Proportion 

x F eature Selectio n v S ignature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE DECOMP (NN,A»UL) 

II. SOURCE (Include coder /per eon with knowledge of program usage): 

N.W. Naugle, Math. Dept. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360 / 65 

IV. PURPOSE: To compute the LU— decomposition of a square nonsingular 
matrix A. Part of package to solve systems of equations and 
compute inverses. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See LU-Decomposition Package, p. 1 of 5 

Description; Decomposes A into LU, where L is lower triangular with ones 
along the diagonal, and U is upper triangular. Uses scaling and 
Gaussian elimination with partial pivoting. 

VI. SPECIAL FEATURES: 

Original matrix A is saved. 

VII. SUBROUTINES REQUIRED: 

SING 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

xPDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : x LACIE x Non-LACIE 

Trainin g x C lassificatio n x E stimation of Proportion 

_x Feature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t y T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE SOLVE (NN,UL,B,X) 

II. SOURCE (include coder /person with knowledge of program usage): 

N.W. Naugle, Math. Dept. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360/65 

IV. PURPOSE: Given the LU-decomposition of A, and a constant vector 
B, the program solves LUX=B for X 

V. ALGORITEM/TECHNIQUE (short description; Include reference if available) 
Reference: See p. 1 

Description: Given a square matrix A with decomposition A=LU, 
and a constant vector B, the solution X of the equation LUX=B 
is obtained. 


VI. SPECIAL FEATURES: can obtain A" 1 by solving LUX.=E. , where E . is a 
unit vector with 1 in i?h slot, i= 1,2,...,NN L 1 1 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x A vailable Not Availabl e W ill be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General put pose^^ Special purpose Program package 

x ppp 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIB 

Trainin g y C lassificatio n v E stimation of Proportion 

x Feature Selectio n y S ignature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE IMPRUV (NN, A,UL,B,X, DIGITS) 

II. SOURCE (include coder /person with knowledge of program usage) : 

N.W. Naugle, Math, Dept. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360/65 

'IV. PURPOSE: To improve the solution X of AX=B 

V. ALGORITHM/TECHNIQUE (short description; include reference If available) 
Reference: See p. 1 

Description: Uses iterative procedure to improve solution X 
(of inverse) of AX=B. The procedure is machine dependent in that 
the tolerance and number of iterations specified depend on 
the word length of the machine- 

VI. SPECIAL FEATURES; 

Double precision accumulation of intermediate sums is required. 

VII. SUBROUTINES REQUIRED: 

SING 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available . 

IX. x S o' 1 ' .ws available off-the-shelf. Approximate Cost: 

I s x I t nut compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : x LACIE x Non-LAC IE 

Trainin g x C lassificatio n y E stimation of Proportion 

x F eature Selectio n x Signature Extension 

Temporal Sampling Systems Developmen t y T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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* * 

I. NAME: SUBROUTINE SING (INHY) A 

II. SOURCE (Include coder /person with knowledge of program usage) : 

N.W. Naugle, Math. Dept. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360/65 

IV. PURPOSE: To print error messages for DECOMP and INPRUV 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See p. 1 

Description: Error messages are printed if original matrix 
in DECOMP is singular or else no convergence is obtained in IMPRUV. 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine 1'nguage conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose x Special purpose Program package 

x ppp H/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Trainin g x C lassificatio n x E stimation of Proportion 

_x Feature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE GINV2M 

II. SOURCE (include coder/person with knowledge of program usage): 

F.M. Speed, Dept. Math. Mississippi State Univ. 

III. PROGRAMMING LANGUAGE: FORTRAN IV, g COMPUTER: 360 

IV. PURPOSE: To compute the generalized inverse of a matrix. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: Adi Ben-Israel and Thomas N. E. Greville, Generalized Inverse — 
Theory and Applications , Wiley Interscience, 1974. 

Description: Uses Gram-Schmidt Orthogonal ization 

VI. SPECIAL FEATURES: Generalized inverse is computed in location 
occupied by original matrix. 


VII. SUBROUTINES REQUIRED: FUNCTION DOT 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: -?_ 

I s x I s not compatible with Image 100 configuration in present form 

Non-FDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION: x LACIE x Non-LACIE 

Trainin g x C lass If icatio n x E stimation of Proportion 

x F eature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: FUNCTION DOT (double precision) 

II. SOURCE (include coder/person with knowledge of program usage) : 
F.M. Speed, Dept. Math. Mississippi State Univ. 

III. PROGRAMMING LANGUAGE: FORTRAN IV, G COMPUTER: 360 

IV. PURPOSE: Inner product of two columns of a matrix. 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 


Description: 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCU MEN TATION, x A vailable Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: __ 

la x is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify _____ _ 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x General purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Trainin g x C lassification x Estimation of Proportion 

x Feature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) — 
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I. NAME: SUBROUTINE CHLSKY (double precision) 

II. SOURCE (include coder/person with knowledge of program usage) : 

Paul R. Hendrick/L. F. Guseman, Jr. Dept. Math. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: XBM. 360/65 

IV. PURPOSE: obtain the modified (sometimes called "alternate form") 

Cholesky decomposition of a real symmetric positive definite 
matrix A, and then its inverse. 

V. ALGORITHM /TECHNIQUE (short description; include reference if available) 


VI. 


VII. 


Reference: Wilkinson, J. H, and Relnsch, C., Handbook for Automatic 
Computation, Vol. II, Linear Algebra, Springer-Verlag, 

New York, 1971, pp. 9-30. 

Description: 

The routine is a translation and slight modification of the 
Algol routine "syminversion" of above reference. The matrix 
A is decomposed into LDL T , with L unit lower tridiagonal and 
D positive diagonal. 


SPECIAL FEATURES: 

The parameter DET input into Chlsky calls for the inverse and/or 
determinant to be computed if DET is non-negative or positive, 
respectively. DET returned as non-positive indicates A is not 
positive definite, either originally or through accumulation of 

roundoff error. 

SUBROUTINES REQUIRED: 


VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x General purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC I E x N on-LACIE 

Trainin g x C lassification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: TSTMIS (Single precision) 

TSTMID (Double precision) 

II. SOURCE ^include coder /person with knowledge of program usage): 

L.F. Guseman,Jr. 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Program for testing matrix Inversion routines for accuracy. 


V. ALGORITHM/TECHNIQUE (short description; Include reference if available) 

Reference: Randall E. Cline, A class of matrices to test Inversion 
procedures, Comm. ACM 7(1964), 724-725. 

Description: Given an nxn matrix T=(t^), where t. .=1 if i#j 


iJ 


t..=d, where d? 4 1, d? 5 -(n-1), the inverse is given by T“ 1 =(1T 1 ..), where 

VI. SPECIAL FEATURES: 

Eigenvalues are given by Xi = (n + d- l) and X-^ =d-l, 2^i^n, and 
can be used to construct matrices with desired condition number. CALL 
Statement must be modified to call routine being tested. 


VII. SUBROUTINES REQUIRED: Inversion routine being tested. 


VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. ^Software available off-the-shelf. Approximate Cost: -?- 

I sx Is not compatible with Image 100 configuration In present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x ppp H/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIEx Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Samplin g S ystems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) . 
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COMPUTER PROGRAM IDENTIFICATION FORM 


1 ot 6 


I. NAME: SUBROUTINE SEIGEN(A,W, AW, C,R,N, DIGIT, NMAX) 

II. SOURCE (include coder/person with knowledge of program usage) : 

W. L. Morris, Math Dept., University of Houston, Houston, Texas 77004 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108/EXEC 8 

IV. PURPOSE: 

Compute eigenvectors and eigenvalues for a real symmetric matrix 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 

Reference: w. L. Morris, Inclusion theorems for a section of a matrix, 
Numer. Math 8(1972), pp. 454-464. 

Description: 

Reduces calculations to that for eigenvalues of sequence of matrices 
of order 2 

VI. SPECIAL FEATURES: 

Requires approximate no. of equivalent decimal 

numbers in floating point representation for the given machine. 

VII. SUBROUTINES REQUIRED: 

SUPSUM RHA 

ORDER RINPRD 

MIND EX 

VIII. DOCUMENTATION: X A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: -?- 

I s X I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

X G eneral purpose Special purpos e x P rogram package 

X PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: X LACIE X Non-LACIE 

Trainin g x C lassificatio n x E stimation of Proportion 

x F eature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t X T echnique Development 

C rop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE MINDEX(R, JM,IM,N) 

II. SOURCE (Include coder /person with knowledge of program usage) : 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: n 0 8 

IV. PURPOSE: Indices of ascending elements of a vector are placed in 
another vector 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 

Description: JM vector will contain the indices of ascending elements 

of R 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: Available X Not Available Will be Available 

IX. X S oftware available off-the-shelf . Approximate Cost: 

I s XI s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

X PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: X LACIE xNon-LACIE 

Trainin g x C lassificatio n x E stimation of Proportion 

x F eature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) , 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: FUNCTION RHA(A,B,C,U,N) 

II. SOURCE (Include coder/person with knowledge of program usage) : 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 1108 

IV. PURPOSE: For a fixed column of A and B matrices, calculate inner 
product A = C * B 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 

Description: u ses rinprd for inner product 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

RINPRD 

VIII. DOCUMENTATION: Available X Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: xLACIE x Non-LACIE 

T rainin g xC lassif icatio n x E stimation of Proportion 
x F eature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t y T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: FUNCTION RINPRD(A,B,N) 

II. SOURCE (include coder/person with knowledge of program usage) : 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: no8 

IV. PURPOSE: Real inner product of vector with another vector 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 

Description: 

A * B = RINPRD 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

SUPSUM 

VIII. DOCUMENTATION: Available x Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non- LAC IE 

Trainin g x C lassificatio n x E stimation of Proportion 

_x Feature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t v T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) _ 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: DOUBLE PRECISION FUNCTION SUPSUM(A, I,N) 

II. SOURCE (Include coder/person with knowledge of program usage) : 

III. PROGRAMMING LANGUAGE: FORTRAN IV, G COMPUTER: 3 60/ 44 

IV. PURPOSE: Vector elements summed 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 

Description: 

VI. SPECIAL FEATURES: Reduced round-off error 

VII. SUBROUTINES REQUIRED: 

ORDER 

VIII. DOCUMENTATION: Available X Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

I s y I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify _ 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

T rainin g K C lassificatio n y E stimation of Proportion 
x F eature Selectio n x S ignature Extension 

Temporal Samplin g S ystems Development v Technique Development 

Crop Assessment Yield Estimation 

Other (specify) ___ 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE ORDER (A, II, 12, N) 

II. SOURCE (include coder/person with knowledge of program usage): 

III. PROGRAMMING LANGUAGE: FORTRAN IV, G COMPUTER: 360/44 

IV. PURPOSE: Orders elements of a vector in increasing absolute value 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 

Description: 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: Available X Not Available Will be Available 

IX. X Software available off-the-shelf. Approximate Cost: 

I s X Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

, Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e X Special purpose Program package 

X PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : X LAC IE XNon-LACIE 

Trainin g x C lassificatio n x E stimation of Proportion 

x F eature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t X T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: 


SUBROUT INE D S VD ( A , MMAX , NMAX , M , N , P , WITHU , WITHV , S , U , V ) 




II. 

III. 


: I 


IV. 


SOURCE (include coder/person with knowledge of program usage) : 

G. H. Golub, Computer Science Dept., Stanford Univ., Stanford, Ca. 94305 
PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 1108, 360/67, CDC6600 

ALGOL 60 

PURPOSE - 

& * Singular value decomposition of rectangular matrix 


j 5 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

i 

Reference: G, H. Golub and C. Reinsch, "Singular value decomposition and 

least squares solutions," Numer. Math 14(1970), pp. 403-421. 

Description: 

^ Decompose a given matrix A into matrices U,S,V such that 

•• A = U * S * V 1 


VI. SPECIAL FEATURES: a may be partitioned such that the submatrix B area 

will have the submatrix U T * B. U and/or V may be 
used in the calculations depending on the flag 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

■ I s X I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

j Some machine language conversion necessary 

In language other than FORTRAN. Specify 

j Software must be designed and developed for Image 100 


X. IMPLEMENTATION LEVEL 

X G eneral purpose Special purpose Program package 

X PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : X LAC IE XNon-LACIE 

Trainin g x C lassificatio n y E stimation of Proportion 

x F eature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t XT echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE MATMUL(A,X,B,I, J,K,NA,NX,NB) 

II. SOURCE (include coder/person with knowledge of program usage): 

David Lowell 

III, PROGRAMMING LANGUAGE: FORTRAN IV, G COMPUTER : 360/44 

IV. PURPOSE. Matrix multiplication 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 
Reference: 

Description; 

A * X = B 

VI. SPECIAL FEATURES: 

Reduces round-off error 
Variable dimensioning 

VII. SUBROUTINES REQUIRED: 

VIII, DOCUMENTATION: Available X Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: ~ 

I s x ls not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language -conversion necessary 

_In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

G eneral purpos e x Special purpose^ Program package 

X PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : x 1ACIE X Non-LAGIE 

Trainin g x C lassificatio n v E stimation of Proportion 

x F eature Selection x Signature Extension 

Temporal Sampling Systems Developmen t y T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE EIGEN 

II. SOURCE (include coder/person with knowledge of program usage): 

Paul R. Hendrick/L. F. Guseman, Jr. Department of Math. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To find all eigenvalues and eigenvectors of a real general 
matrix. 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: Smith, B. T. , et al., Matrix Eigensystem Routines — EISPACK 
Guide, Springer-Verlag, Lecture Notes in Computer Science, 
vol. 6, New York, 1973. 

Description: This is merely a driver routine and calls in order the 

routines which balance the matrix, transform it to upper 
Hessenberg form and accumulate the elementary similarity 
transformations, compute the eigenvalues and eigenvectors of 
the balanced matrix, and then backtransform these to the 

VI. SPECIAL FEATURES: eigenvectors of the original matrix. 

If the parameter IERR is returned as nonzero, the eigenvalue with its 
value as index was not computed in 30 iterations. Eigenvalues indexed 
from IERR + 1 to N should be correct. Eigenvectors have not been computed 
in this case. 

VII. SUBROUTINES REQUIRED: 

BALANC, BALBAK, ELMHES, ELTRAN, HQR2 ,' UNPK 

VIII. DOCUMENTATION: x Available Not Available Mill be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x l a not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


B.l-20 


COMPUTER PROGRAM IDENTIFICATION FORM 2 of 7 

I. NAME: SUBROUTINE BALANC 

II. SOURCE (include coder/ person, with knowledge of program usage) : Burton S. Garbov, 
Argonne National Laboratory, Applied Mathematics Division, 9700 S Cass Ave. , 

III. r^^^^IING 1 ?Sg!&AGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: This subroutine balances a real matrix and isolates eigenvalues 
whenever possible. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: Smith, B.T., et al.. Matrix Eigensystem Routines-EISPACK Guide, 
Springer-Verlag, Lecture Notes in Computer Science, New York, 1973. 

Description: 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration In present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software 'must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Tr aining C lassification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Samplin g S ystems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME; SUBROUTINE BALBAK 

II. SOURCE (include. coder/person with knowledge of program usage) : Burton S. Garbow, 
Argonne National Laboratory, Applied Mathematics Division, 9700 S. Cass Ave. , 

III. PR0g§&eS 1 M^A^S; 39 FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: This subroutine forms the eigenvectors of a real general matrix by 
back transforming those of the corresponding balanced matrix determined by 
BALANC. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: Smith, B.T. , et al.. Matrix Eigensystem Routines-EISPACK Guide, 
Springer-Verlag, Lecture Notes in Computer Science, vol. 6, New York, 1973. 

Description: 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x i s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

_x General purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

__Temporal Sampling Systems Development Technique Development 

C rop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 4 of 7 

I. NAME: SUBROUTINE ELMHES 

II. SOURCE (Include coder /person with knowledge of program usage) : Burton S. Gar bow, 
Argonne National Laboratory, Applied Mathematics Division, 9700 S. Cass Ave, 

III. ^IS?SXiai5i 1 ISS^9A&S1 39 FORTRAN IV COMPUTER: IBM 360/65 

IV. PURP0SE:Given a real general matrix, this subroutine reduces a submatrix 

to upper Hessenberg form by stabilized elementary similarity transformations. 

V. ALGQRITHM/TECHNIQUE (short description; include reference if available) 

Reference :Smith, B.T., et al., Matrix Eigensystem Routines-EISPACK Guide, 
Springer-Verlag, Lecture Notes in Computer Science, vol. 6, New York, 1973. 

Description: 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION; x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

_x_PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) . 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE ELTRAN 


II. 

III. 


SOURCE (include coder /person with knowledge of program usage) : Burton S. Garbow 
Ar S onne N filinois L 60439 t0ry ’ Applied Mathematics Division, 9700 S. Cass Ave. , ’ 
PR0GRA^liiiG IANGUAGE: 3 9 FORTRAN IV COMPUTER: IBM 360/65 


IV. PURPOSE: This subroutine accumulates the stabilized elementary similarity 
transformations used in the reduction of a real general matrix to upper 
Hessenherg form by ELMHES. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: Smith, B.T., e£ al., Matrix Eigensystem Routlnes-EISPACK Guide, 
Springer-Verlag, Lecture Notes in Computer Science, New York, 1973. 

Description: 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: ; 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpos e Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

C rop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROCRAM IDENTIFICATION FORM 
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I. 

II. 

III. 

IV. 


NAME: SUBROUTINE HQR2 


SOURCE (Include coder/person with knowledge of program usage) : Burton S. Garbow 
N fii?^L L |g?§| tor y» Applied Mathematics Division, 9700 S. Cass Ave., 


FORTRAN IV 


COMPUTER: IBM 360/65 


PURPOSE: This subroutine finds the eigenvalues and eigenvectors of a real 
upper Hessenberg matrix by the QR method. 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference :smith, B.T., et al. , Matrix Eigensystem Routines-EISPACK Guide, 
Springer-Verlag, Lecture Notes in Computer Science, vol.6, New York, 1973. 

Description: 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x A vailable _Not Availabl e W ill be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose _Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE UNPK 

II. SOURCE (include coder/person with knowledge of program usage): 

Paul R. Hendrick/L. F. Guseman, Jr. Dept. Math. TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To break up the complex eigenvector matrix returned by 

eigenvector routines into its real and imaginary parts. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: Smith, B. T., et al.. Matrix Elgensystem Routines — EISPACK 
Guide, Sprlnger-Verlag, Lecture Notes in Computer Science, 
vol. 6, New York, 1973, pp. 73-75. 

Description: The eigenvectors are originally packed in one matrix, one 
real eigenvector per column corresponding to a real 
eigenvalue, or the real part and imaginary part in two 
adjacent columns corresponding to a conjugate pair of 
eigenvalues . 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x l a not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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GENERAL MATHEMATICAL LIBRARY 
2. OPTIMIZATION PACKAGE 


COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE DFMFP 

II. SOURCE (include coder/person with knowledge of program usage): 

IBM 3 60 /Scientific Subroutine Package-L. F. Gusernan, Jr . 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360/65 

IV. PURPOSE: To find a local minimum of a function of several variables 

by the method of Fletcher and Powell. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: R. Fletcher and M.J.D. Powell, A rapidly convergent descent 

method for minimization, Comput. J. 6(1963), 163-168 
Description: 

Uses deflected-gradient method. Converges in a finite number 
of steps if the objective function is a quadratic form. 

VI. SPECIAL FEATURES: 

Requires a user-supplied subroutine which provides the values 
of the objective function and the gradient for each argument vector. 

VII. SUBROUTINES REQUIRED: 

FUNCT - User supplied subroutine as described above. 

VIII. DOCUMENTATION; x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I b not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must he designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION ; x LACIE x Non-LACIE 

Training Classificatio n x E stimation of Proportion 

x F eature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t y„ T echnique Development 

C rop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE FRKWLF 


A II. 


III. 


SOURCE (include coder/person with knowledge of program usage) : 
B. C. Peters, Jr. / L. F. Guseman, Jr. j Bruce Marion 
PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: i B M 360/65 


IV. 


- V. 


PURPOSE: To solve the quadratic programming problem : minimize | | Ax-y | [ 

subject to r , 

f 1 x 4 = ^ , x. > 0, 1 < i < n 
i=l 1 1 ~ ~ ~ 

ALGORITHM/TECHNIQUE (short description; include reference if available) 


Reference: M. Frank and P. Wolfe, An Algorithm for Quadratic 
Programming, Naval Research Logistics Quarterly, Vol. 3, 1956. 

Description: A modification of the Frank-Wolfe algorithm 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED; 

VIII. DOCUMENTATION: Available x Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: ~2z. 

I s xl s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non- LAC IE 

Training Classificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) ; 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I . NAME: SUBROUTINE QM431 

1 1 . SOURCE (include coder /person with knowledge of program usage) : 

Guseman, Jr. / Bruce Marlon 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Minimization of a quadratic form subject to linear 
inequality and nonnegativity constraints on the variables 
(Quadratic programming problem) 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: A. Ravindran, Algorithm 431-A computer routine for 
quadratic and linear programming problems, Comm. ACM. 15 (1972), 
818-820. 

Description: Solves the quadratic program: minimize z=c^x+x^Qx 
subject to Ax£fb, x>0 using Lemke's complimentary pivot algorithm 
Can also be used to solve linear programming problems. 

VI. SPECIAL FEATURES: Original matrix Q is lost. 


VII. SUBROUTINES REQUIRED: INITIA, NEWBAS, SORT, PIVOT 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x Is not__ compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE INITIA 

II. SOURCE (include coder/person with knowledge of program usage): 

L . 1 ; . Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To find the initial almost complementary solution by 
adding an artificial variable 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 
Reference: See QM431 

Description: 

VI. SPECIAL FEATURES: Part of QM431 package 


VII. SUBROUTINES REQUIRED: 




VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classificatio n x E stimation of Proportion 

Feature Selection ^Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) . 
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X, NAME: SUBROUTINE N^t^BAS 

II. SOURCE (include coder/person with knowledge of program usage): 
L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To find the new basis column to enter in terms of the 
current basis in performing a simplex-type algorithm 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See QM431 


Description; 


VI. SPECIAL FEATURES: Part of QM431 package 


VII. SUBROUTINES REQUIRED: PPRINT 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify _______________ 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Training Classificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE PPRINT 

II. SOURCE (include coder/person with knowledge of program usage): 

L.F, Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To place the current solution to complementary problem 
into the array ALPHA 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See QM431 

Description: 

VI. SPECIAL FEATURES: Part of QM431 Package 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

_x_PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Training Classificatio n x E stimation of Proportion 

Feature Selection_ Signature Extension 

emporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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5 oi 6 


I. NAME: SUBROUTINE SORT 

II. SOURCE (include coder/person with knowledge o£ program usage): 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To find the pivot row for next iteration by the use of 
(simplex-type) minimum ratio rule. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See QM431 

Description: 

VI. SPECIAL FEATURES: Part of QM431 package 

VII. SUBROUTINES REQUIRED : PPRINT 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Sof tware must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Trainin g C lassificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE PIVOT 

II. SOURCE (include coder/person with knowledge of program usage): 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To perform the pivot operation by updating the inverse 
of the basis and Q vector. 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 
Reference: See QM431 

Description: 

VI. SPECIAL FEATURES: Part of QM431 package 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : x LACIEx Non-LACIE 

Training Classificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 


I. NAME: NCIMEN 

II. SOURCE (include coder /person with knowledge of program usage) : 

Henry Decell, University of Houston, Math Dept. 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: UNIVAC 1108 

IV. PURPOSE: T° determine the simultaneous solution of an overdetermined 

system of nonlinear equations. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: An algorithm for solving overdetermined system of non 

linear equations, John Engvall MSC-IN-67-ED-49 (Januarv, 
1967). 

Description; 

Problem is set up in least squares fashion and solution 
is computed by use of the Newton-Raphson iterative technique. 

VI. SPECIAL FEATURES: There must be at least as many equations as unknowns. 


VII. SUBROUTINES REQUIRED: Subroutine MOPS must be supplied by user to 

evaluate partial derivatives in addition to 
an initial guess at solution. 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. y S oftware available off-the-shelf. Approximate Cost; _n_ 

I s X I s not compatible with Image 100 configuration in present form 

X Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e X Special purpose Program package 

X PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-IACIE 

T raining Classificatio n X E stimation of Proportion 

X F eature Selectio n X S ignature Extension 

Temporal Sampling Systems Developmen t X T echnique Development 

C rop Assessment Yield Estimation 

Other (specify) 
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UUMTUXER rROUKAM IDENTIFICATION FUKM 


x or y 


I. NAME: Basic Spline Package - Single Variable 

II. SOURCE (Include coder /person with knowledge of program usage) : 

Philip W. Smith, Dept, of Math. TAMU 
III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To handle all computations involved in evaluating B-spline 
functions and their derivatives, including conversion to piecewise 
polynomial representation. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Ref erence: Package for calculating with B-splines, Carl de Boor, 
the University of Wisconsin-Madison, Mathematics Research Center, 

MRC Technical Summary Report #1333, October 1973. 

Description: See above reference. 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: BSPLDR, BSPLEV, BSPLPP, BSPLVN, BSPLVD, BVALUE, 
INTERV, PPVALU 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: ~ ~ 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpos e x P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

F eature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessmen t Y ield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE BSPLDR 

II. SOURCE (include coder/person with knowledge of program usage) : 

Philip W. Smith 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Constructs divided difference table for B-spline coefficients 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See reference for Basic Spline Package-Single Variable 

Description: Equation 2.(2), above reference, is employed. 

VI. SPECIAL FEATURES: Cases which would require division by zero 

because of coincidence of knots are not evaluated: their values are 
bypassed in the calling routines. 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selectio n S ignature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 


B.3-3 



COMPUTER PROGRAM IDENTIFICATION FORM 


3 of 9 


I. NAME: SUBROUTINE BSPLEV 

II. SOURCE (Include coder/person with knowledge of program usage) : 

Philip W. Smith 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Calculates value of the B-spline and its derivatives 
(up to a specified order) at a particular x value. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference; See references for Basic Spline Package-Single Variable 

Description; Equation 2.(1), above reference, is used. 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: INTERV , BSPLVN 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : x LAC IE x Non- LAC IE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE BSPLPP 

II. SOURCE (include coder/person with knowledge of program usage) : 

Philip W. Smith 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Converts B-spline representation to piecewise 
polynomial representation. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See reference for Basic Spline Package-Single Variable 

Descriptions The necessary parameters are calculated by BSPLDR and 
BSPLEV . 


VI. SPECIAL FEATURES: Requires temporary storage "scratch" array allocated 
in calling program. 


VII. SUBROUTINES REQUIRED: BSPLDR, BSPLEV 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. Software available off-the-shelf. Approxi ma te Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purposex Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Training ^Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

C rop Assessmen t . Y ield Estimation 
Other (specify) - 
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I. NAME: SUBROUTINE BSPLVD 

II. SOURCE (Include coder /person with knowledge of program usage) : 

Philip W. Smith 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Calculates value and derivatives (up to specified 
order) of all B-splines which do not vanish at the x value. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See reference for Basic Spline Package-Single Variable 

Description: Repeated calls are made to routine BSPLVN and an array 
is built with the different order derivatives at the different pieces 
of the domain. 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: BSPLVN 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e r X Special purpose Program package 

x P DP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : x LACIE x Non-LACXE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE BSPLVN 

II. SOURCE (include coder /person with knowledge o£ program usage): 

Philip W. Smith 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Calculates the value of all possibly non— zero B-splines 
of the particular order (either input or previously determined) 
at the x value. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See reference for Basic Spline Package-Single Variable 

Description: See p.20 above reference to algorithm which this routine 
incorporates. 


VI. SPECIAL FEATURES: An argument parameter switch can be set to use 
previously computed values where this can save time. 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 


B.3-7 



COMPUTER PROGRAM IDENTIFICATION FORM 


7 of 9 


I. NAME: FUNCTION BVALUE 

II. SOURCE (include coder/person with knowledge of program usage) : 

Philip W. Smith 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Calculates value at x of the specified-order derivative 
from the B-representa tion. 

V. ALGORITHM/TECHNIQUE (short description; Include reference, if available) 
Reference: See reference for Basic Spline Package-Single Variable 

Description: See section 2, above reference. 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: INTERV 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: _ 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x_ Special purpose Program package 

x PDP H/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

T emporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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x. NAME: subroutine interv 

II. SOURCE (include coder /person with knowledge of program usage): 

Philip W. Smith 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. • PURPOSE: Determines largest subscript for a sequence such that 
the Indexed term is less than a specified value x. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See reference for Basic Spline Package-Single Variable 

Description: See p.24, above reference. 


VI. SPECIAL FEATURES: An attempt is made to use the previous call's 
results to minimize search time. 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is_ x I s not compatible with Image 100 configuration in present form 

Non-PDF 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: FUNCTION PPVALU 

II, SOURCE (include coder/perscra with knowledge of program usage) : 

Philip W. Smith 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV, PURPOSE: Calculates value at x of a specified order derivative 
from the piecewise polynomial representation. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See reference for Basic Spline Package-Single Variable 

Descriptions Equations 1.(1), above reference, are used. 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: INTERV 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Co3t: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x ppp 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIEx Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development fechnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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B. GENERAL MATHEMATICAL LIBRARY 
4. TRANSFORM PACKAGE 

A fairly general program which provides the required transform 
capabilities is available from PAR, Inc. 
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I. NAME: (MSFAR) Multispectral Fourier Analysis Routines 

II. SOURCE (Include coder /person with knowledge of program usage): 

Jack D . Bryant 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: C DC 6400 

IV. PURPOSE: To extract spatial features from multispectral scanner data 
sets. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: J. A. Smith, G. Hartinger and T. Johnson, "Investigation of 
Preprocessing Transformations", final report, vol.III, NAS9-12972, 

Dec. 1974. 

Description: See reference above, volumes I and II, also. 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

x S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

y G eneral purpose Special purpos e v P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE xNon-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessmen t Y ield Estimation 

_y_0ther (specify) Spatial RecngniMnn Pattern 
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I. NAME: MERGE 

II* SOURCE (include coder /person with knowledge of program usage): 

Jack D. Bryant 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTES: C0C 6400 

IV. PURPOSE: To read up to 12 data files and merge selected spectral 
channels from each file onto one file. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See MSFAR references above. 

Description: Data handling 


VI. SPECIAL FEATURES: Examples: the data set may consist of a large 

number of channels (e.g.24 channels) or may be a large file (e.g. 
an ERTS image) . 


VII. SUBROUTINES REQUIRED: MOVER* GETD, PUTD* ’ SKIPF* 
*Written is COMPASS SETC* GNUM* 

(CDC 6400 Assembly Language) CHARS* ZFET* 


VIII o DOCUMENTATION: g Available Not Available Hill be Available 


IX. x S oftware available off-the-shelf. Approximate Cost: 

la Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

x S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 


X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

* 

XI. APPLICATION; LACXE x Non-LACIE 

Training Classification Estimation of Proportion 

JFeature Selection Signature Extension 

^Temporal Sampling Systems Development Technique Development 

Crop Assessmen t ' Y ield Estimation 

x O ther (specify) Spatial Pattern Recognition 
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I. NAME; CHNL 

II. SOURCE (Include coder /person with knowledge of program usage): 

Jack D . Bryant 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER; CDC 6400 

IV. PURPOSE: To build a RECOG formatted data file from values read from 
punched cards for use by the Fourier transform program. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See MSFAR reference, above 

Description; Data handling 

VI. SPECIAL FEATURES: 




VII. SUBROUTINES REQUIRED: PUTD*, SETC* 

*Written in COMPASS 

(CDC 6400 Assembly Language) 

VIII. DOCUMENTATION: x A vailable Not Availabl e W ill be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

x S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed far Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selectio n S ignature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

x O ther (specify) Spatial Pattern Recognition 
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I. NAME: FFT2D 

II. SOURCE (Include coder /person with knowledge of program usage): 

Jack D. Bryant 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: CDC 6400 

IV. PURPOSE: FFT2D obtains fast Fourier transforms of the input data streams 
employing dynamic dimensioning for major arrays. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See MSFAR reference, above 

Description: Driver program for FOURT. (Fast Fourier Transform Subroutine) 


VI. SPECIAL FEATURES: Two modes of operation are provided: only one 

aperture within the flight line can be used or the window can be forced 
to move through the flight line. 

VII. SUBROUTINES REQUIRED: DOIT, CORE*, SETC*, CALCP, FOURT, SKIPF* 

RELOC, GETO, PUTD* 

*Written in COMPASS (CDC 6400 Assembly Language) 

VIII. DOCUMENTATION: x Available Not Available Mill be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

x S ome machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

x O ther (specify) Spatial Pattern Recognition 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE F0UR2 

II. SOURCE (Include coder/person with knowledge of program usage) : 

Jack. D. Bryant 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: IBM 360 

IV. PURPOSE: To obtain the Cooley-Tukey Fast Fourier Transform 

V. ALGORITHM/ TECHNIQUE (short description; Include reference if available) 
Reference: IEEE Audio Transactions, June 1967, special issue on FFT. 

Descriotion: Obtains two-dimensional FFT when the number of data points is 
2 n * 2 * 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 


VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration In present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 


X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 



APPLICATION : LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

T emporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 


x Other (specify) Spatial Pattern Recognition 
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I. NAME: SUBROUTINE F0UR1 

II. SOURCE (include coder/person with knowledge of program usage): 

Jack D . Bryant 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: IBM 360 

IV. PURPOSE: To obtain the Cooley-Tukey Fast Fourier Transform 

V. ALGORITBM/TECHNIQUE (short description; include reference if available) 
Reference: IEEE Audio Transactions, June 1967, special issue on FFT 

Description: Obtains one-dimensional FFT when the number of data points 
is a power of two. 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 


VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approx ima te Cost: _____ 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 


X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

T emporal Sampling Systems Development Technique Development 

v Crop Assessment Yield Estimation 

x O ther (specify) Spatial_ Pattern, Recognition 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE FOURT 

II. SOURCE (include coder/person with knowledge of program usage): 

Jack D. Bryant 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: CDC 6400 

IV. PURPOSE: To obtain the Cooley- Tukey Fast Fourier Transform. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: IEEE Audio Transactions, June 1967, special issue on FFT. 

Description: This is the fastest and most versatile version of the 
FFT known. 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x ppp H/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

'" v Crop Assessment Yield Estimation 

_x_Other (specify) Spatial Pattern Rarngm.t-frm 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: TRANSF 

II. SOURCE (Include coder /person with knowledge of program usage): 

Jack D. Bryant 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: CDC 6400 

IV. PURPOSE: Driver program for user— supplied program TRANS which samples 
data generated by FFT2D and writes a RECOG formatted tape or file. 

This may then be analyzed using RECOG or CLUSTD. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: See MSFAR reference 


Description : Data handling 


VI. SPECIAL FEATURES: 


VII. 


SUBROUTINES REQUIRED: GETD , PUTD* , TRANS (user supplied) ,C0RE* , SETC* , 
PTAPE 


* Written in COMPASS (CDC 6400 Assembly Language) 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDF 11/45 FORTRAN 

x S ome machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIEx Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t T echnique Development 

C rop Assessment Yield Estimation 

x O ther (specify) Spatial Pattern Recognition 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: NORMLX 

II. SOURCE (include coder/person with knowledge of program usage) : 

Jack D. Bryant 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: CDC 6400 

IV. PURPOSE: Normalizes and rewrites in RECOG format ouput of the circular 
TRANS subroutine. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See MSFAR reference 

Description : 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: READTP, GETD, PUTD* 

*WRITTEN IN COMPASS (CDC 6400 Assembly Language) 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

x S ome machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpos e P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : LAC IE x Non- LAC IE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

jCrop Assessment Yield Estimation 

x O ther (specify) Spatial Pattern Recognition 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: WSORT 

II. SOURCE (include coder/person with knowledge of program usage): 

Jack D. Bryant 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: CDC 6400 

IV. PURPOSE: Reads and performs straight or nomalized sort of spatial 
features obtained by wedge sampling, and produces a new data file, 
in RECOG format, of sorted spatial features. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: See MSFAR reference. 

Description: 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: SORT, NSORT, RDTAPE, WRTAPE, GETO, PUTD* 

*WRITTEN IN COMPASS (CDC 6400 Assembly Language) 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

x Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

x Other (specify) Spatial Pattern Recognition 
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c. GENERAL STATISTICAL LIBRARY 
1. STATISTICAL SUMMARY 


C.l-1 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: SUBROUTINE GROPl (X , N , NG , XMIN , XMAX , F ) 

II. SOURCE (include coder /person with knowledge of program usage): 

Theory & Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: 1108 

IV. PURPOSE: 

Group data into specified number of equally spaced groups 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Description: Calculate frequencies of each point in each group 

VI. SPECIAL FEATURES: XMIN = lower end of the first group 

XMAX - upper end of the last group 
NG = number of groups 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose^ Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) __ 
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COMPUTER PROGRAM IDENTIFICATION FORM 1 of 1 

I. NAME: SUBROUTINE MEDIAN(X,N,XMED) 

II. SOURCE (include coder /person with knowledge of program usage) : 

UNIVAC STAT PACK 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

IV. PURPOSE: Calculate median of a sequence of data points 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: JSC STAT CAT 

UNIVAC STAT PACK 

Description: Computes rank of median and observations 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is r.jt compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI, APPLICATION: LAC IE Non- LACIE 

Training Classification Estimation of Proportion 

Filature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) _ 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE MODE (X,N,XM0D ? XB,XM2,XM3,XM4,G1,G2,SKEW,K) 

II. SOURCE (include coder/person with knowledge of program usage) : 

JSC Theory of Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

IV. PURPOSE: Calculate estimate of mode of a distribution from sample 

of data 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: jsc STAT CAT 

Description: 


VI. SPECIAL FEATURES: Finds central moments, correct 3rd and 4th central 

moments, calculate coefficient of skewness and excess. 
Distributions must be unimodal and belong to class 
of Pearson system distributions (gamma, normal) 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessmen t Y ield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I* NAME 1 SUBROUTINE RANGE (X,N,R) 

II. SOURCE (include coder /person with knowledge of program usage): 

UNIVAC STAT PACK 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 110 8 

IV. PURPOSE: Calculate range of sequence of numbers 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Description: Range = X - x , 

max min 

VI. SPECIAL FEATURES: 


VII. SUM. 1ES REQUIRED: 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE Non-LACIE 

Trainin g C lassification Estimation of Proportion 

Feature Selection Signature Extension 

T emporal Sampling Systems Development Technique Development 

C rop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE MNDEV (X,N, IND,D) 

II. SOURCE (Include coder/person with knowledge of program usage) : corrected from 
JSC Theory and Analysis Office, Houston STAT PACK 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

IV. PURPOSE: Calculate mean deviation of an array of data 


V, ALGORITHM/TECHNIQUE (short description; include reference if available) 


Reference: JSC STAT CAT 

UNIVAC STAT PACK 


VI. 


Description: 

deviation = ^ J |x.-u| u = median or mean 

i=l 1 X. => i 1 " ■ observation 

x 

N = no. of observation 

SPECIAL FEATURES: 

Option to find deviation from median or arithmetic mean 


VII. SUBROUTINES REQUIRED: MEDIAN (if deviation to be found from median) 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE JSfon-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling_ Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE SDEV(X,N s IND,XB,S) 

II. SOURCE (include coder /person with knowledge of program usage) : 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN ¥ COMPUTER: 1108 

IV. PURPOSE: Calculate standard deviation of sequence of data points 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Descript ion: 

„ _ v6c 2 -{SX) 2 /N . „ _ v/IX 2 -(£X) 2 /N 

s or s s 

VI. SPECIAL FEATURES: Option if estimate to be used in unbiased or maximum 

likelihood 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION; X Available Not Available Will be Available 

IX. X Software available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : LAC IE Non-LACIE 

Trainin g C lassification Estimation of Proportion • 

Feature Selectio n S ignature Extension 

T emporal Samplin g S ystems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: SUBROUTINE ORDER (X,N) 

II. SOURCE (include coder /person with knowledge of program usage) 
JSC Theory and Analysis Office, Houston 
III, PROGRAMMING LANGUAGE: FORTRAN V COMPUTER. 1108 


IV. PURPOSE: Arrange an array in ascending order 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC .STAT CAT 

Description: 


VI. SPECIAL FEATURES: Ordered data returned in X array locations 


VII. SUBROUTINES REQUIRED: 


VIII. 


IX, 


DOCUMENTATION :_X_Available Not Available Will be Avail 

X Software available off-the-shelf. Approximate Cost. , ■ - 

IT_Ib not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify — 

Software must be designed and developed for Image 100 


X. IMPLEMENTATION LEVEL 

General purpose_Special purpoee_Program package 
POP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling_Systems Development_Technique Developmen 

Crop Assessment Yield Estimation 

Other (specify) " 
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COMPUTER PROGRAM IDENTIFICATION FORM 
THIS FORM PREPARED BY LEC 


1 of 1 


I. NAME: PLOT PAC 

II. SOURCE (include coder /person with knowledge of program usage): 

Ted Kell 

III, PROGRAMMING LANGUAGE: FORTRAN/MACRO COMPUTER: 11/ 45 

IV. PURPOSE: To provide the user with a variety of data plotting options 
(histogram /scatter gram/ graphs/etc. ) 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 
Reference: 

Description; 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION : Available Not Available Will be Available 

IX. Software available off-the-shelf. Approximate Cost: _ 

Is Is aot w __coapatihle with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE y Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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THIS FORM PREPARED BY LEG 

I. NAME; Binary Search/N-Dimensional Histogram Acquisition 

II. SOURCE (include coder /person with knowledge of program usage); 

General Electric/Jim Brierly 

III. PROGRAMMING LANGUAGE; FORTRAN COMPUTER; Image-100/PDP-ll 

IV. PURPOSE; To reduce the time neccessary to acquire N-Dimensional 
histograms. 

V. ALGORITHM/ TECHNIQUE (short description; Include reference if available) 
Reference: See next page 

Description; See next page 


VI. SPECIAL FEATURES; See next page 


VII. SUBROUTINES REQUIRED :None 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: S6.300 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpos e x P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION ; xLACIE xNon-LACIE 

x T rainin g x C lassification Estimation of Proportion 

x F eature Selectio n x S ignature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 


C.l-io 
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1. 2,2 BINARY SEARCH/N-DBIENSIONAL HISTOGRAM ACQUISITION 

The ti^r.e necessary to acquire N-dimcnsioual histograms is normally significantly 
reduced using this program. A binary search technique is employed to (potentially) 
reduce the number of cell interrogations required for signature generation by testing 
large volumes of feature space rather than each individual cell. This is accomplished 
by successively bisecting each feature space axis, each time testing each produced 
volume for non-zero pixel counts. Any volume with a zero pixel count is subsequently 
ingored, and in most cases will enable a substantial reduction in the number of single 
cells otherwise requiring consideration. Any \'olume with a non-zero count will be 
further subdivided until the required resolution is reached. The pixel counts for the 
smallest cells are recorded as is done in the N-D Histogram. 

Use of this program generally results in ordors-of-magnitude reductions in histogram 
acquistion time, resulting in N-dimensional training being more attractive as an 
analysis tool. In addition, the operation of post-training functions such as semi-super- 
vised cluster analysis ( described in paragraph 1. 4. 2 ) * becomes more efficient. 


< 


< 


C'.l-ll 


c. GENERAL STATISTICAL LIBRAE? 

2. DENSITIES & DISTRIBUTIONS PACKAGE 


C.2-1 


COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: FUNCTION PHINV (P) 

XI. SOURCE (include coder /person with knowledge of program usage): 

JSC THeory and Analysis Office, Houston 
III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: no8 

IVi PURPOSE: Inverse of normal distribution function 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Description; Given mean 0 and variance 1, using Newton Raphson 
iterations 

VI. SPECIAL FEATURES: Checks probability P for validity and prints error 

message if P is invalid. 


VII. SUBROUTINES REQUIRED: 

PHI 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) ... __ 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. 


II. 

III. 


NAMEl DOUBLE PRECISION FUNCTION DPHI(X) 

SOURCE (toc^fuife ^co<?er99erson. with knowledge of program usage) : 
S.W. Kalin g , JSC/LEC, Houston 

PROGRAMMING LANGUAGE-; FORTRAN V COMPUTER: u 


IV. PURPOSE: Normal probability integral for real X 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: jsc STAT CAT 


Description: 

1_ 

X&i 



exp (“t 2 /2)dt 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION : x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Samplin g S ystems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 
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I. NAME: REAL FUNCTION FISH(F,N1, N2) 


II. SOURCE (include coder /person with knowledge of program usage): 


corrected from UNIVAC STAT-FACK: JSC Theo 

III. PROGRAMMING LANGUAGE: FORTRAN ' 


sis Office. Houston. 
108 


IV. PURPOSE: Approximate Fisher's f-distribution 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

JSC STAT CAT 

Reference :Qeorge Neel, Larry Whitehead, Robert Bottenburg, "Calculation 
of the percentage points of the F-distribution," Technical Documentary Report 
PFL-TDR-64-16, 6570th Personal Research Lab., Aerospace Medical Div., 
Description: Air Force Systems Command, Lackland 

AFB , Texas 


VI. SPECIAL FEATURES: N1 and N2 > 1 


VII. SUBROUTINES REQUIRED: PH£ 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


C.2-4 



COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: DOUBLE PRECISION FUNCTION DERF(X) 

, , , FUNCTION ERF (X) , . 

II. SOURCE (include coder /person with knowledge of program usage): 

S.W. Kabng, JSC/LEC, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: n08 

IV. PURPOSE: Error function 

V. ALGORITBM/TECHNIQUE (short description} Include reference if available) 
Reference: jsc STAT CAT 

Description: 

f exp ( - t 2 ) dt 
Vf T ) 

J o 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf . Approximate Cost: 

Is Is not compatible *rith Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software muBt be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

C rop Assessment Yield Estimation 

Other (specify) 


C. 2-5 



COMPUTER PROGRAM IDENTIFICATION FORM 


3 of 1 


I. 

II. 

III. 

IV. 

V. ALGORITHM/ TECHNIQUE (short description; include reference If available) 
Reference: jsc STAT CAT 

Description: l-ERF(K) 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: v Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost; 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

. Temporal Samplin g S ystems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) . 


DOUBLE PRECISION FUNCTION DERFC(X) 

, . . , .FUNCTION ERFC(X), , 

SOURCE (include coder/peraon with knowledge of program usage) : 

S.W. KAHNG, JSC/LEC, Houston 

PROGRAMMING LANGUAGE : FoRTRAN v COMPUTER: U08 

PURPOSE: Complimentary 
Error function 


C.2-6 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. 

II. 

III. 

IV. 

V. 


VI. 

VII. 

VIII. 

IX. 


X. 

XI. 


NA 1 ® 1 REAL FUNCTION CHI2 (X,V) 

SOURCE (include coder /person with knowledge of program usage): 
JSC Theory and Analysis Office, Houston 
PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 110 8 

PURPOSE: Approximate distribution 


ALGORITHM/TECHNIQUE (short description; include reference If available) 
Reference: JSC STAT CAT 

Description: Given a point and degrees of freedom, approximate value 

of the distribution at the point 

Method used depends on values of the point and the degree of freedom. 

SPECIAL FEATURES: Provision for overflow 

X > 0.0 

uses CBRT function 


SUBROUTINES REQUIRED: ZIP (calls FACTOR) 

OVERFL 

PHI 

DOCUMENTATION: X Available Not Available Will be Available 

X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify , 

Software must be designed and developed for Image 100 

IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

APPLICATION : LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessmen t Y ield Estimation 

Other (specify) 


C.2-7 




COMPUTER PROGRAM IDENTIFICATION FORM 1 of 1 


I. NAME; FUNCTION EXPAN (N,T) 

II. SOURCE (Include coder /person with knowledge of program usage) : 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: n 0 8 

IV. PURPOSE: Approximate student’s t-distribution 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: jsc STAT CAT 

Description: Using series expansion method, given degrees of freedom 


VI. SPECIAL FEATURES: Values of N decide which equations to use. 


VII. SUBROUTINES REQUIRED: 

I 

VIII. DOCUMENTATION : XAva liable Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


C. 2-8 



COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME; REAL FUNCTION POIS(N,XMU) 

II. SOURCE (include coder /person with knowledge of program usage); 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: u 0 8 

IV. PURPOSE: Value of Poisson distribution function 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

N i 

Description: POIS(N,XMU) = e -XMU j ; XMU = mean of distribution 

i=0 i ! 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X Software available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) . 


C.2-9 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: FUNCTION FISHIN (ALPHA, NX, N2) 

II, SOURCE (Include coder/person with knowledge of program usage) : 

Corrected from UNIVAC STAT PACK 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: uo8 

IV. PURPOSE: Inverse Fisher T s distribution 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: jsc STAT CAT 

Description: ref. Abramowitz and Stegun, Handbook of Mathematical 
Function , p.947. National Bureau of Standards. 

VI. SPECIAL FEATURES: 0<ALPHA<1, ALPHA is the confidence coefficient. 


VII. SUBROUTINES REQUIRED: FISH PHINV 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

r 

IX. X S oftware available off-the-shelf . Approximate Cost: 

la Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify __________________ 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : LACIE Non-LACIE 

T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

T emporal Sampling Systems Development Technique Development 

C rop Assessment Yield Estimation 

Other (specify) ___ 


C, 2-10 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: real FUNCTION BIN(N,P,M) 

II. SOURCE (include coder/person with knowledge of program usage): 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

IV. PURPOSE: Value of binomial distribution function 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: jsc STAT CAT 

Descriptio: : 


VI. SPECIAL FEATURES: For values of M > 1000, PHI is called in the 

calculations 


VII. SUBROUTINES REQUIRED: 

PHI 

ZOT 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


C . 2-11 



COMPUTER PROGRAM IDENTIFICATION FORM 


1 of X 


I. NAME: FUNCTION HYTRIC(K,N,N1,NR) 

II. SOURCE (include coder /person with knowledge of program usage): 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: -nos 

IV. PURPOSE: Hypergeometric distribution 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Description: Calculate probability that a random variable with a 

hyper geometric distribution is less than or equal to 
a given value 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

FCTELG 

FACTOR 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Samplin g S ystems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) __ 


C. 2-12 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: REAL FUNCTION EETINC(X,A, B) 

II. SOURCE (Include coder /person with knowledge of program usage) : corrected from 

UNIVAC STAT PACK 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

IV. PURPOSE: INCOMPLETE BETA DISTRIBUTION 


V. ALGORITHM/TECHNIQUE (short description; Include reference if available) 

Reference: JSC STAT CAT 

UNIVAC STAT PACK 

Description! Compute probability that a random variable, distributed 
as the incomplete Beta function with given values of 
parameters, is less than or equal to a given value. 


VI. SPECIAL FEATURES: The Modified Jacobi Quadrature formula is used in 

computing integrals. Uses CBRT function 


VII. SUBROUTINES REQUIRED: 

PHI 

CHI2 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 


C , 2-13 


COMPUTER PROGRAM IDENTIFICATION FORM 


I of 1 


I. NAME: real FUNCTION GAMIN (X, A) 

II. SOURCE (include coder /person with knowledge of program usage): 
UNIVAC STAT PACK 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

IV. PURPOSE: Incomplete gamma function 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: JSC STAT CAT 

UNIVAC STAT PACK 


Description: 

P(A,X) = / e -t t A_1 dt with P N (A,X) 

where N is such that ^N-1^ A ’^ ^ 1rt -6 

VI. SPECIAL FEATURES: - 10 

A and X > 0 


rA e "X y X n 
* n i 0 r (A+n+1) + 


VII. 'SUBROUTINES REQUIRED: 


VIII. DOCUMENTATION: X A vailable Not Available Will be Available 

IX. xS oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Samplin g S ystems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


C . 2-14 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: REAL FUNCTION STUDIN(N,P) 

II. SOURCE (include coder /person with knowledge of program UBage) : 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: no8 

IV. PURPOSE: Inverse student’s t- value 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: jsc STAT CAT 

Description: Given degrees of freedom and significance level, using 
Newton Raphson iterations 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

FACTOR PHINV 
EXPAN 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available Gf f-the-shelf . Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) ' 


C. 2-15 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of L 


I. NAME: REAL FUNCTION ZOT(N) 

II. SOURCE (include coder/person with knowledge of program usage) : 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

IV. PURPOSE: ZOT « ALOG(Nl) 

V. ALGORITHM/TECHNIQUE (short description.; include reference if available) 
Reference: JSC STAT CAT 

Description: Calls FACTOR with a = N + 1 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

FACTOR 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE Non-LACIE 

Training Classificatio n E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


C. 2-16 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. HAKE: FUNCTION RANDOM 

II. SOURCE (include coder/person with knowledge of program usage): 

Richard Rosencranz, JSC Engineering Applications, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: n 0 8 

IV. PURPOSE: Generate random numbers from uniformly distributed set 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Description: 


VI. SPECIAL FEATURES: Option to return number from uniformly distributed 

set or number that is a function of input 
Used mainly as a random starting point for cyclical set of numbers 
in Z0R function. 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. _X Software available off-the-shelf. Approximate Cost: 

I s I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X . IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Samplin g S ystems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: REAL FUNCTION ZIP (A) 

II. SOURCE (include coder /person with knowledge of program usage): 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: H08 

IV. PURPOSE: Compute ln^ of the gamma function 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Descriptions Calls FACTOR with a = A 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

FACTOR 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify ___ 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


C.2-18 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: REAL FUNCTION ZARF(Z) 

XI. SOURCE (include coder/person with knowledge of program usage): 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

IV. PURPOSE: Compute complete gamma function T(Z) 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Description: Calls FACTOR with a = Z + 1.0 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

FACTOR 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Samplin g S ystems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


C.2-19 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: FUNCTION KPSNSM (AS) 

FUNCTION KPSNLG (AL) 

II. SOURCE (include coder/person with knowledge of program usage): 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE : FORTRAN V COMPUTER : n08 

IV. PURPOSE: random numbers from a Poisson distribution. 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

k+1 

Description: KPSNSM: jt V. < 

1=1 x 
k+1 

KPSNLG: E In <— y 

i=l 

VI. SPECIAL FEATURES: use KPSNLG if AL is greater than 88,028. If Z0R 
is initialized with a large odd number in the calling program, the 
sequence of random numbers will be repeated. 


VII. SUBROUTINES REQUIRED: Z0R 

VIII, DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

C rop Assessment Yield Estimation 

Other (specify) , 


C. 2-20 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: SUBROUTINE BIN2 (N, I, BINC9),D7) 

II. SOURCE (Include coder/person with knowledge of program usage): 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: uo8 

IV. PURPOSE: Compute binomial coefficient. 

V. ALGORITHM/TECHNIQUE (short description; Include reference if available) 
Reference: jsc STAT CAT 

Description: Comppte N things taken 1 at -a time. 

(*) D N! 

\I/ I! (N-I) ! 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED* Values of I and N determine the calculations done. 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose_ Special purpose Program package 

PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection _Signature Extension 

Temporal Samplin g S ystems Development Technique Development 

C rop Assessment Yield Estimation 

Other (specify) 


C. 2-21 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: SUBROUTINE OVERFL(I) 

II. SOURCE (include eoder/person with knowledge of program usage) : 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: H08 

IV. PURPOSE: Checks values for overflow condition 

V, ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 

Description: Checks for overflow condition and return a one if 
overflow occurred 

VI. SPECIAL FEATURES: OVERFL is part of 1108 system 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: Available X Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) __________ 


C. 2-22 


C. GENERAL STATISTICAL LIBRARY 
3. REGRESSION PACKAGE 


GENERAL STATISTICAL LIBRARY 


4. STATISTICAL TESTS 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of I 


I. NAME; SUBROUTINE TTEST (A, N , M, NN , MM, X, AL , IND , S , US , AF, AT) 

It. SOURCE (Include coder /person with knowledge of program usage); 
JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE! FORTRAN V COMPUTER: no 8 

IV. PURPOSE; Ho tel ling's T 2 test on means. 


V. ALGORITHM/TECHNIQUE (short description; Include reference if available) 


Reference ; jsc STAT CAT 


Description; Test at a desired confidence level the null hypothesis 
that the mean vector of a set of observations equals a given vector. 


VI. SPECIAL FEATURES; Variable dimensioning of arrays. Acceptance or 

rejection message is printed and option is available to print output. 


VII. SUBROUTINES REQUIRED: C0VARI GINV2M FISHIN 

VIII. DOCUMENTATION; v Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

T emporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) „ 


C.4-2 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: SUBROUTINE C0VARI (A s N,M,MM,NN } S,X) 

II. SOURCE (Include coder/pereon with knowledge of program usage): 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE * F0RTRAN v COMPUTER: 1108 

IV. PURPOSE: Compute covariances between observed variables. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Description: 

( i \ N M M r 

S(1 5 J) = ^N-lj SEE A(K,I) - X(I)| j A(K } J) - X(J)1 

K-l 1=1 >1 L J L J 

VI. SPECIAL FEATURES: 

A is a matrix where column i contains the ovservations for variable 
i . N and M -5 2 

l 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpos e P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

T emporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


C.4-3 




COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: SUBROUTINE ANOVl(Y,N,J , I,R,S,SS ,NT>F, SM, F,PF, $K,IND, YH) 

II. SOURCE (Include coder/person with knowledge of program usage); corrected from 

JSC Theory and Analysis Office, Houston UNIVAC STAT PACK 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 110g 

IV. PURPOSE: One way analysis of variance 


V. ALGORITHM/TECHNIQUE (short description; include reference If available) 
Reference: JSC STAT CAT 


Description: One way analysis of variance with unequal number of 
replications in each cell; computes all sums of squares, degrees of 
freedom, mean squares, F-levels and probabilities that F-levels are 
exceeded. 

VI. SPECIAL FEATURES: Control may be returned to a labelled statement in 

calling program in event of overflow 
Option flag specifies lines to be printed. 


VII. SUBROUTINES REQUIRED: FISH 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. _X Software available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image .100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


C.4-4 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: SUBROUTINE ANOVA2(Y,N,NI,NJ,NK,B,R,SS,NDF.,SM,F,PF, $K, JFR,YH) 

II. SOURCE (include coder/person with knowledge of program usage) : corrected from 

JSC Theory and Analysis Office, Houston UNIVAC SI AT PACK 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: n 08 

IV. PURPOSE: Two way analysis of variance 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Description: 2 way analysis of variance with an equal number of 

replications in each cell; computes all sums of squares, 
degrees of freedom, mean squares, F-levels and probabilities 
that F-levels are exceeded. 

VI. SPECIAL FEATURES: If overflow occurs, control may be returned to some 

labelled statement in the calling program 
Option flag determines type of output desired. 

VII. SUBROUTINES REQUIRED: FISH 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. XS of tware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) ___ . 
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COMPUTER PROGRAM IDENTIFICATION FORM 


I of I 


I. NAME: SUBROUTINE RMNOV(Y,N,NB,NT,NR,B,R,SS,NDF,SM,F,PF, $ K) 

II. SOURCE (include coder/person with knowledge of program usage) : 
JSC Theory Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

IV, PURPOSE: Analysis of variance on a randomized block design 


V. ALGORITHM/TECHNIQUE (short description; include reference If available) 
Reference: JSC STAT CAT 


Description: Analysis of variance on a randomized block design with an 
equal number of replications in each experimental unit; 
computes all sums of squares, degrees of freedom, mean 
squares, F-levels and probabilities that the F levels are 
exceeded. 

VI. SPECIAL FEATURES: If overflow occurs, control will be returned to a 

labelled statement in the calling program 


VII. SUBROUTINES REQUIRED: 

BTALOG 

PROBLY 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

S ome machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: FUNCTION ZOR(N) 

II. SOURCE (include coder/peraon with knowledge of program usage): 

Richard Rosencranz, JSC Engineering Applications, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 11G 8 

IV. PURPOSE: Random number generator 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Description: Uses uniform distribution set between 0 and 1 or normal 
distribution with mean 0 and variance 1 

VI. SPECIAL FEATURES: Flag to choose distribution set 

Calling program should declare ZOR to be an ABNORMAL FUNCTION 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X Software available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

C rop Assessment Yield Estimation 

Other (specify) __ 
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COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAMES SUBROUTINE KOLSMR (X,N,F,KW,KN s D) 

II. SOURCE (include coder /person with knowledge of program usage) : 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

IV. PURPOSE: Kolmogorov-Smirnov test of goodness of fit 


V. ALGORITHM/TECHNIQUE (short description; Include reference if available) 
Reference: JSC STAT CAT 

Description: Test goodness of fit of data to a specified theoretical 
cumulative distributive function 


VI. SPECIAL FEATURES: Cumulative distributive function may be applied to 

raw or normalized data; prints output 
Print option specifies lines to be printed out. 

Function option specifies type of function F is to be used. 

VII. SUBROUTINES REQUIRED: 

ORDER 

DOCUMENTATION: X Available Not Available Will be Available 

X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image IQO 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


VIII. 

IX. 
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COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 1 


I. NAME: SUBROUTINE FACTOR (ALPHA, EXPAN, YEXPAN) 

II. SOURCE (include coder /person with knowledge of program usage): 

JSC Theory and Analysis Office, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: nog 

IV. PURPOSE: Approximate complete gamma (factorial) function 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: JSC STAT CAT 

Description: Evaluate gamma function and In of the complete gamma 
function 

VI. SPECIAL FEATURES: Value of alpha determines equations used. 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: X A vailable Not Available Will be Available 

IX. X Software available off-the-shelf. Approximate Cost: 

la Is no t c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC1E Nou-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

__ Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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D. PATTERN RECOGNITION LIBRARY 


1 


CLUSTERING PACKAGE 


COMPUTER PROGRAM IDENTIFICATION FORM 
THIS FORM PREPARED BY LEC 


EPPLER 
1 of 2 


I. NAME: Semi-supervised Cluster Analysis 

II. SOURCE (include coder /person with knowledge of program usage) : 

General Electric/ Jim Brier ly 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: IMAGE-100/PDP-11 

IV. PURPOSE; To determine the feature space location of homogenous clusters 
from the results of N-dimensional histogram analysis. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See next page 

Description: See next page 


VI. SPECIAL FEATURES: See next page 


VII. SUBROUTINES REQUIRED: None 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: $6,049 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpos e x P rogram package 

PD? 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION; x LACIE xNon-LACIE 

x T rainin g x C lassification Estimation of Proportion 

x F eature Selectio n x S ignature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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1. 4.2 SEMI-SUPERVISED C LU STER ANALYSIS 

The function of this program is to determine. « passible the • 

space location of homogeneous cell groupings (c ’- 1S eF S)' _ , results of a previous 

Iona! ■histogram analy sis, in a sen.*- supervised »”“ F - «"■? ?! O v ! ,1, ■ ”o 
N-d'i«casionai mMK step, SMS program MWM *® eatesorwc *c> -O 
a user-specified number .«£ clusters. A migrating weans teclwnjue t* »sw for part. 
tlZZ feature space with the weans M*My set to avtatravy swr«K nte. The 
program Aon iterates M a succession «mm ***** to deternnne the o,*n.um 
asstatment of each cell to one of *e nuwber of clusters specdied by the user, ho 
stop Ju pMcrrm. at any point awl accept the resuilts to-date, or may elect 
to let the prolan. run to coneinafon. In either ease, .the program displays Hie results 
• on (he l«sptay.Ter»M. » is possible that an acceptable sulu ton may 

not ho reached for (hr given ^Mt «r any »«h« «««*« of speetted clusters, the 
program \vll!i tciiu-Uvatc aftev Iviity 


cn 
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COMPUTER PROGRAM IDENTIFICATION FORM 
THIS FORM PREPARED BY LEC 


1 of 3 


I. NAME: DRIVER PROGRAM (ISOCLS) 

II. SOURCE (include coder /person with knowledge of program usage): 

R. MINTER/LEC also E. KAN/LEC * 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: UNIVAC 1108 EXEC 2 

SLEUTH 

IV. PURPOSE: Driver program to perform iterative clustering 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: "Computer Program Documentation, ISOCLS, Program C094" 

R. Minter, LEC, Report #CPD202, Oct. 1972. "User Documentation 
EOD-LARSYS, Program Q6I9" R. Minter, LEC, Report //CPD406, 
Description: LEC-3984, Nov. 1974. 

Iteratively group data into homogenous nominal-sized clusters, 
the nominal size is to be preset by users. 

VI. SPECIAL FEATURES: See references 


VII. SUBROUTINES REQUIRED: TDATE, SETUP, TAPERD, ISODAT, COVARR, PCHSTA, CHAIN, 

PRINT, DASTAP (See next page) 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. _x Software available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image ICO configuration in present form 

Non-PDF 11/45 F0RTR'\N 

Some machine language conversion necessary 

_In language other ban FORTRAN. Specify 

Software must be designed and developed for Ima'ge 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

x T raining Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 


s-ra 


General f low of the? Program ~£ so <=-^ <, 


♦ 



Obtains Reads card 
current input and 
date makes drum 

assignments 
FLTNUM 
NXTCHR 
NUMBER 
DRMAVL 1 
DRMASG 2 


Reads tape 
input and 
stores on 
drum. 
UNPACK 
FRMTl 
FRMT2 
NTRAN 3 4 


Performs 

clustering 

algorithm. 

PRINT 

CLDIST 

NTRAN 3 


Calculates 
and prints 
covariance 
matrices 
NTRAN 3 


Punches 
means and 
covariance 
matrix 


Performs Prints Creates 
chaining final* 3 tape for 
algorithm results PMIS-DAS 
NTRAN 3 PACK 
SETMRG 4 NTRAN 3 


1 System routine which returns the number of words available on FH432 and FH1782 drums. 

2 System routine which assigns specific word length to drum unit. 

3 System binary I/O routine. 

4 System routine which sets the margin at the top and bottom of the page. 
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COMPUTER PROGRAM IDENTIFICATION FORM 


3 of 3 


I. NAME: SUBROUTINES TDATE, SETUP, etc. (See first page of ISOCLS) 

II* SOURCE (include coder/person with knowledge of program usage) : 

III. PROGRAMMING LANGUAGE: COMPUTER: 

IV. PURPOSE: 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference : 

Description: 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: Available Not Available Will be Available 

IX. Software available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Samplin g S ystems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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D. PATTERN RECOGNITION LIBRARY 


2. CLASSIFICATION PACKAGE 


Maximum likelihood classification algorithms are being 
implemented on the Image 100 System by LEC. 
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COMPUTER PROGRAM IDENTIFICATION FORM EPPLER 

THIS FORM PREPARED BY LEC 1 of 2 

I. NAME: ELLTAB 

II. SOURCE (include coder /person with knowledge of program usage): 

NASA/MTF Clay Jones 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: UNI VAC 1108 

IV. PURPOSE: Classify multispectral scanner data. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See next page 

Description: See next page 

VI. SPECIAL FEATURES: See next page 

VII. SUBROUTINES REQUIRED: None 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 0 

I s x I s not compatible with Image 100 configuration in present form 

x Non-PDP 11/45 FORTRAN 
Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Trainin g x C lassification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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AN IMPROVED VERSXOJ^QF^T?fE„TABLK LCQK^UP 
ALQ ORlTHM^EOa ^FATTER N RECOGNITION * 

W* G. Eppler 

Lockheed Electronics Company* 

Houston Texas 

ABSTRACT 

The Table Look-Up Approach to pattern recognition has been 
used for 3 years at several research centers In a variety of 
applications. A new version has been developed which is faster, 
requires significantly leas core memory, and retains full pre- 
cision of the input data* The hew version can be used on low* 
cost minicomputers having 32 K words (16 bits each) of core mem* 
ory and fixed-point arithmetic; ho special*purpose hardware is 
required* An initial FORTRAN version of this system can clas- 
sify an ERTS computer-compatible tape into 24 classes In loss 
than 15 minutes* 


iHRLEMENTATION OF AN ADVANCED fABlJl LQOK*UR J^ 
FOB LARGE AREA LAND-UEH C L ASS1F 1 CATlgN 
Olay Jones 


Johnson Space Center 
Earth Resources Laboratory 
National Aeronautics and Space Administration 
Bay St* ‘Louis , Mississippi 39520 


ABSTRACT 

Software employing Ropier* s improved table 
look-up approach to pattern recoct nit ion has been 
developed, and results from this software are 
presented. The look-up table for each class is 
a computer representation of a hypcrollipsoid in 
four dimensional space. During imp! enen tat ion of 
the software Eppler ' s look-up procedure was modi- 
fied to include multiple ranges in order to accom- 
modate hollow regions in the ell ipso ids. In a 
typical ERTS classification run less than 6000 
3 G-bit computer words were required to store tables 
for 24 classes. Classification results from the 
improved table look— up arc identical with those 
produced by the conventional method, i . e . , by cal^ 
culation of the maximum likelihood decision rule 
at the moment of c lass if i cat ion • With the tew 
look-uo approach an entire LRTS MSS frame can be 
classified into 24 classes in 1.3 hours, compared 
to 22.5 hours required by the conventional method. 
The new software is coded completely in FORTRAN to 
facilitate transfer to other digital computers. 
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COMPUTER PROGRAM IDENTIFICATION FORM bPPLER 

THIS FORM PREPARED BY LEC 

I. NAME : Maximum Likelihood Resolution 

II. SOURCE (include coder /person with knowledge of program usage): 

General Elec trie /Jim Brier ly 

III, PROGRAMMING LANGUAGE: FORTRAN COMPUTER: IMAGE-100/PDP-11 

IV. PURPOSE: To resolve conflicts between overlapping spectral signatures. 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 
Reference: See next page 

Description: See next page 

VI. SPECIAL FEATURES: See next page 


VII. SUBROUTINES REQUIRED: None 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: $6,286 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpos e x P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION: xLACIEx Non-LACIE 

Trainin g x C lassification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) __ 
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1.3.1 MAXIMUM LIKELIHOOD RESOLUTION 


This program permits conflicts between overlapping spectral signatures to be resolved 
by assigning spectral cells to the class which has highest probability , where probability 
Is based on the frequency of occurrence of a given multtspoctrul pixel value normalized 
With respect to die traiiflng area. 


\ - 

The user first specifies 'two signature files (results of N- Dimensional Histogram 
Acquisition) presumably generated from the same image data* The option of simultaneous 


theme modification Is available and if desired, the one or two corresponding themes 
tracks arc also specified. The image data must also be stored in the Image Memory 
Unit. 


The two signature files are compared in feature space, and any areas of overlap are 
examined. The overlapping cells arc assigned to one or the other signature on the 
basis of maximum likelihood. That is, for any feature space cell, the normalized 
number of pixels having that cell's coordinates is computed for each signature, and 
the largest one is retained in its file while the other file's cell is deleted; thus, 
assuring that each cell is assigned to only one signature file. The themes, if desired, 
are also modified accordingly by deleting the pixels corresponding to the deleted cells. 


D. PATTERN RECOGNITION LIBRARY 


3. FEATURE SELECTION PACKAGE 


COMPUTER PROGRAM IDENTIFICATION FORM 


1 of 14 


I* NAME: PROGRAM UH Linear Feature Selection 

II. SOURCE (include coder /person with knowledge of program usage): 

J. A. Quirein, LEC-JSC Houston; H. P. Decell/R. Teun, Math. Dept., UH 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: 1108/EXEC 8 

360/67 

IV. PURPOSE: Linear feature selection to minimize loss of separabilicy 
among classes. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: H.P. Decell & J.A. Quirein, "An Iterative Approach to the Feature 
Selection," Univ. of Houston Report 26, March 1973, NAS 9-12777. Kenneth 
Baker, "Users Guide to the Univ. of Houston Linear Feature Selection Program," 
Description: NASA, 1973. Rebecca Teun, "Addendum to Users Guide," UH Report 33, 

July, 1974 

Linear transformation B of measurement x such that the NAS 9-12777 

average transformed divergence is maximized for all 
interclass pairs. 

VI. SPECIAL FEATURES: Reduces dimensions of data to be classified from 

n to k. Uses -David on Iterator to find best B. 


VII. SUBROUTINES REQUIRED: See tree on p. la. 

DAVIDN SHUTLE MPROD EIGEN TRANDV 
CUBIC CHLSKY XPROD BHATCH 

PINT MVPR HALM DMVPR 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X Software available off-the-shelf. Approximate Cost: 

Is Is no t X c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose ^ Special purpose X Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION: XLACIE XNon-LACIE 

Training Classification Estimation of Proportion 

X F eature Selection Signature Extension 

Temporal Samplin g S ystems Development Technique Development 

C rop Assessment Yield Estimation 

Other (specify) _ 
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COMPUTER PROGRAM IDENTIFICATION FORM 2 of 14 


I. NAME: SUBROUTINE DAVIDN 

II. SOURCE (include coder/person with knowledge of program usage): 
Ivan Johnson, JSC/MPAD, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: 1108/EXEC B 

i 60/67 

IV. PURPOSE: Davidon-Fletcher-Powell iterator 


V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 

Reference: Ivan Johnson, "Impulsive Orbit Transfer Optimization by an 
Accelerated Gradient Method, "Journal of Spacecraft and Rockets, Vol. 

6, No. 5, May 1969. 

Description :Davidon- Fletcher-Powell optimization method to find best 
k linear combination B matrix. 


VI. SPECIAL FEATURES: Optimal B is derived iteratively from an initial 
guess B and gradient. 


VII. SUBROUTINES REQUIRED -.CUBIC, FINT, SHUTLE 

VIII. DOCUMENTATION: Available x Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

la Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling _3ystems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 


i ol L4 


I. NAME: SUBROUTINE FINT (FNT, IPART) 

II. SOURCE (include coder/person with knowledge of program usage): 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

360/67 

IV. PURPOSE: Compute numerical partial derivatives. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See DAVIDN 

Description: 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: SHUTLE 

VIII. DOCUMENTATION: Available x Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t y c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: y LACIE xNon-LACIE 

Training Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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COMPUTER PROGRAM IDENTIFICATION FORM 


4 of 14 


I. NAME: SUBROUTINE CUBIC (XX, YY,XMIN2) 

II. SOURCE (include coder/person with knowledge of program usage): 

Ivan Johnson, JSC/MPAD, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

360/67 

IV. PURPOSE: Cubic fit technique. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See DAVl'DN 

Description: 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: Available x Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configu- i - n in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non- LAC IE 

Training Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE SHUTLE (FNT , IP ART) 

II. SOURCE (include coder/person with knowledge of program usage): 

J. Quirein, LEC/JSC; H.P. Decell, R. Teun, Math. Dept. UH 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: nos 

360/67 

IV. PURPOSE: Calculate interclass divergences and probability of misclassifi- 
cation bound. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: J.A. Quirein, "Iterative Approach to Feature Selection Problem," 
UH Report 26, March 1973, NAS9-12777 

Description: See above report for equations used. 

VI. SPECIAL FEATURES : Option to find Bhattacharyya bound and transformed 
average divergence. 


VII. SUBROUTINES REQUIRED: MPROD, CHLSKY, HALM, BHATCH, TRANDV 

VIII. DOCUMENTATION: Available x Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

G eneral purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION ;x LACIEx Non-LACIE 

Training Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE CHLSKY (XTEMP, XUP, CHI, K) 

II. SOURCE (include coder/person with knowledge of program usage): 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: 1108 

360/67 

IV. PURPOSE: Invert a positive definite symmetric matrix. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 

Description: XTEMP is factored into 2 upper triangular matrices that are 
inverted. 


VI. SPECIAL FEATURES: Option to check inverse by multiplying the inverse by the 
original matrix and printing the result. 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: Available x Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Trainin g x C lassification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE MVPR(XTEMP, VECT, XVECT, KDM, KN) 

II. SOURCE (include coder/person with knowledge of program usage) : 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

360/67 

IV. PURPOSE: Post-multiply mxn matrix by n*l vector. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 

Description jXTEMP x VECT = XVECT 


VI. SPECIAL FEATURES: Single precision. Dimensions less than or equal to 12. 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: Available x Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x 1 s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-IACIE 

Trainin g C lassification Estimation of Proportion 

Feature Selection Signature Extension 

T emporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) . 
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I. 

II. 

III. 

IV. 

V. 


VI. 

VII. 

VIII. 

IX. 


X. 

XI. 


NAME: SUBROUTINE MPROD (XX,YY,ZZ,K1,K2,K3) 

SOURCE (include coder /person with knowledge of program usage) : 

PROGRAMMING LANGUAGE ; FORTRAN V COMPUTER: 1108 

36U/67 


PURPOSE: Matrix multiplication. 


ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference : 

Description: XXxYY = ZZ 


SPECIAL FEATURES : Double Precision: Dimensions must be less than OR 
equal to 12. 


SUBROUTINES REQUIRED: 

DOCUMENTATION: A vailabl e x N ot Available Will be Available 

x S of tware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

APPLICATION : x LAC IE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE XPR0D(XX,YY,ZZ,K1 S K2,K3) 

II. SOURCE (Include coder/person with knowledge of program usage) : 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

360/67 

IV. PURPOSE: Single precision matrix multiplication, 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 
Reference: 

Description: XX*YY=ZZ 


VI. SPECIAL FEATURES: Dimensions must be less than or equal to 12. 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: Available x Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify __ 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP H/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) _____ ___ 
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I. NAME: SUBROUTINE HALM(N, IFLAG, ERFT) 

II. SOURCE (include coder/person with knowledge of program usage) : 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

360/67 

IV. PURPOSE: Pairwise bounds on probability of misclassif ication. 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 
Reference: See references for UH feature selection. 

Description: 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: EIGEN, XPROD, MVPR 

VIII. DOCUMENTATION: Available x Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify _ 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAClEx Non-LACIE 

Trainin g x C lassification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) __ 
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I. NAME: SUBROUTINE EIGEN(AA,N,MV,A,E,R) 

II. SOURCE (include coder /person with knowledge of program usage): 

Jane Montgomery, TRW Systems 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

360/67 

IV. PURPOSE: 

eigenvalues and eigenvectors of a real symmetric matrix 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: Von Neumann, "Mathematical Methods for Digital Computers," 
Ed. A. Ralston and H. S. Wilf, John Wiley and Sons, New 
York, chapter 7, 1962. 

Description: 

Jacobi method of diagonalization 

VI. SPECIAL FEATURES: 

option to compute eigenvalues only 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: Available XNot Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

I s X I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

X G eneral purpose Special purpose Program package 

X PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: X LAC IE xNon-LACIE 

Training Classification Estimation of Proportion 

F eature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE BHATCH(F.NT, IPART) 

II. SOURCE (include coder /person with knowledge of program usage): 

J. A. Quirein, LEC/JSC, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: nos 

360/67 

IV. PURPOSE: 

compute minimum average Bhattacharyya bound on the probability of 
misclassif ication 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: UH feature selection 

Description: 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

MPROD 

CHLSKY 

DMVPR 

VIII. DOCUMENTATION: Available XNot Available Mill be Available 

IX. X S oftware available off-the-shelf. Apprpximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose X Special purpose Program package 

X PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: XLACIE xNon-LACIE 

Training Classification Estimation of Proportion 

X F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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i. 

ii. 

in. 

IV. 

V. 


VI. 

VII. 

VIII. 

IX. 


X. 

XI. 


NAME: SUBROUTINE DMVPR(XTEMP, VECT, XVECT, KDM,KN) 

SOURCE (Include coder/person with knowledge of program usage) : 


PROGRAMMING LANGUAGE*. FORTRAN V 
PURPOSE: 

multiply matrix by vector 


COMPUTER: 110g 

360/67 


ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 

UH feature selection 


Description: 

XTEMP * VECT = XVECT 


SPECIAL FEATURES: 

double precision 


SUBROUTINES REQUIRED: 

DOCUMENTATION: Available x Not Available Will be Available 

X S oftware available off-the-shelf. Approximate Cost: 

I s X I s no t c ompatible with Image 100 configuration in present form 
Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

IMPLEMENTATION LEVEL 

X G eneral purpose Special purpose Program package 

X PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

APPLICATION : X LAC IE XNon-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE TRANDV 

II. SOURCE (include coder/person with knowledge of program usage): 

J. A. Quirein, LEC/JSC, Houston 

III. PROGRAMMING LANGUAGE: FORTRAN V COMPUTER: 1108 

360/67 

IV. PURPOSE: 

minimum average transformed divergence 

V. ALGORITHM/ TECHNIQUE (short description; include reference If available) 
Reference: UH feature selection 

Description: 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

MPROD 

CHLSKY 

DMVPR 

VIII. DOCUMENTATION: Available X Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is no t Xc ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

X PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

X F eature Selection Signature Extension 

Temporal Samplin g S ystems Deve.1 opmen t T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: LFSPMC : Linear Feature Selection Probability of Misclassif ication 

II. SOURCE (include coder /person with knowledge of program usage): 

L.F. Guseman,Jr. / Bruce D, Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360/65 


IV. PURPOSE: To compute a single linear combination of features which 
minimizes the probability of misclassif ication 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: L.F. Guseman,Jr. B. Charles Peters, Jr., and Homer F. Walker, 

On minimizing the probability of misclassif ication for linear feature selection, 
Ann. Statist. 3(1975), 1-8 

Description: See program package documentation LFSPMC: Linear Feature Selection 
Program Using Probability of Misclassification, Contract NAS-9-13894-1S, 

Report #3, Dept, of Math. Texas A&M University 


VI. SPECIAL FEATURES: 

COMPLETE PROGRAM PACKAGE: INDIVIDUAL SUBROUTINES AND 
DRIVER ATTACHED. 


VII. SUBROUTINES REQUIRED: 

See subroutine tree on la of 18 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration In present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpos e x P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: DRIVER PROGRAM (LFSPMC) 

II. SOURCE (include coder/person with knowledge of program usage): 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360/65 

IV. PURPOSE: Driver program to set common dimensions throughout 
LFSPMC Package 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See previous LFSPMC 

Description: Sets value MX for maximum numbers of classes and value 
NX for maximum number of features. Provides maximum string lengths. 


VI. SPECIAL FEATURES: MX and NX must not be exceeded during any run for 
a given compilation. 


VII. SUBROUTINES REQUIRED: PRDIM 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Soma machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : x LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

C rop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE PRDIM 

II. SOURCE (include coder/ person with knowledge of program usage): 

L.F. Guseman, Jr. / Bruce Marion. 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To provide dimensioning information to the subroutines of 
LFSPMC, effectively simulating object-t ime dimensioning, 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Ref erence : LFSPMC 

Description: Starting addresses of subset arrays contained within single, 
large, previously-defined one-dimensional arrays are computed. 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

BMAT 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE BMAT 

II. SOURCE (include coder/person with knowledge of program usage): 

L. F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Main Driving Routine for LFSPMC. It performs all input and 
output of class statistics, executes user specified program 
options. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: LFSPMC 

Description: 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

BVECT, CLSNEW, ERRMAT, RANCLN 

VIII. DOCUMENTATION: y Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is not x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpose x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : x LAC IE x Non- LAC IE 

Training Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Samplin g S ystems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE BVECT 

II. SOURCE (include coder/pereon with knowledge of program usage) : 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Compute the vector B which minimizes the probability 
of misclassification function 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 
Reference: LFSPMC 

Description: BVECT is a driver routine. It inputs iteration 
parameters, either Inputs a starting value for the B vector or 
calls for its computation (BCOMP) , calls for iteration on B 
for function minimization computation (DFMFP), and calls for 
computation and output of the final (lower case) function value. (FUNCT. ) 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

BCOMP, DFMFP, FUNCT 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non--PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

T raining Classification Estimation of Proportion 

x F eature Selectio n S ignature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) , 
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I. NAME: SUBROUTINE BCOMP 

II. SOURCE (Include coder/person with knowledge of program usage) : 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER:IBM 360/65 

IV. PURPOSE: Compute the initial B vector for input to the minimization 
iteration routine. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: LFSPMC 

Description: Determines iterative soution of a special case of the 
M-class problem (equal covariances) and uses this as starting vector 


VI. SPECIAL FEATURES: 

No computation is performed in the special case of two classes 
with equal a priori probabilities; it merely calls routine BC2CP. 


VII. SUBROUTINES REQUIRED: 

BC2CP, EIGEN 

VIII. DOCUMENTATION: x Available Not Available Mill be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpos e P rogram package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

x F eature Selectio n S ignature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other . (specify) 
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I. NAME: SUBROUTINE BC2CP 

II. SOURCE (include coder /person with knowledge of program usage): 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Compute initial vector Bo far input to the minimization 

iteration routine - 2 class case with equal a priori probabilities 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Ref erence : LFSPMC 


Description: 

Bo=(p!-p 2 ) T (E 1 +E 2 )~ 1 

VI. SPECIAL FEATURES: 

Called only in the special case of two classes with equal a 
priori probabilities 


VII. SUBROUTINES REQUIRED: 


VIII. 


IX. 


X. 


XI. 


I 


DOCUMENTATION: x A vailable Not Available Will be Available 

x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

APPLICATION : x LAC IE x Non-LACIE 

Training Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) — — — 
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I. NAME: SUBROUTINE EIGFN 

II. SOURCE (include coder /person with knowledge of program usage) : 

L.F. Guseman, Jr. / ‘Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Compute eigenvalues and eigenvectors of a real 
symmetric matrix. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: SYSTEM/360 SSP and 

"Mathematical Methods for Digital Computers," Vol. 1, Edited by 
A. Ralston and H.S. Wilf, John Wiley and Sons, N.Y., 1962, Chapt. 7. 
Description: 

Implements the Jacobi diagonalization method as adapted by 
Von Neumann 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration In present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

X PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE DFMFP 

II. SOURCE (include coder/person with knowledge of program usage) : 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: To find a local minimum of a function nf several variables. 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 

Reference: R. Fletcher and M.J.P. Powell, A rapidly convergent 
descent method for minimization, Comp. J. 6(1963), 163-168. 

Description: 

Uses a deflected gradient procedure to determine a local 
minimum of a function of several variables. Requires analytical 
expressions for gradient. 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

FUNCT - A user supplied routine, which computes the objective 
function and gradient. 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x General purpose Special purpose Program package 

x PDF 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION: x LACIE x Non-LACIE 

Training Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE FUNCT 

II. SOURCE (include coder /person with knowledge of program usage): 

III, PRbGlR^li® ai lM^U^E ; Bruce ^ion COMPUTER: 

FORTRAN IV IBM 360/65 

IV. PURPOSE: Compute the probability of misclassif ication and its 
gradient with respect to a vector 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: LFSPMC 

Description: Uses analytical variational equations and closed 
form expressions for the probability of misclassif ication to compute 
functional values and gradient vector for optimization package 

VI. SPECIAL FEATURES: Has variable dimensioning 


VII. SUBROUTINES REQUIRED: XNDF, DPHI, DERF 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non- LAC IE 

Training _C las sif ication Estimation of Proportion 

x Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME; FUNCTION XNDF 

II. SOURCE (Include coder/person with knowledge of program usage): 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE; FORTRAN IV COMPUTER; IBM 360/65 

IV. PURPOSE; Evaluate one-dimensional normal density function with 
mean XMU and variance SIG at a point A. 

V. ALGORITHM/TECHNIQUE (short description; Include reference if available) 

Reference: LFSPMC:Linear Feature Selection Program Using the 
Probability of Misclassif ication 

Description: XNDF(A,XMU,SIG)=plJT 1 ^ 2 (SIG)“ 1 ^ 2 EXP (-4 


VI. SPECIAL FEATURES: Sets XNDF(A s XMU,SIG)=0 if argument in EXP is 
less than -175 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Training Classification Estimation of Proportion 

x F eature Selectio n S ignature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE DERF 

II. SOURCE (include coder /person with knowledge of program usage): 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: ibm 360/65 

IV. PURPOSE: To evaluate the error function of a variable X. 

V. ALGORITHM/TECHNIQUE (short description* include reference if available) 
Reference: 

X 

Description: DERF(X)= (2tf _l ^ J EXP( ) dt 

-a 

VI. SPECIAL FEATURES: Has additional entry point DPHI to compute 
DPHI (X)*,5 + .3 (DERF(X/V^2~)) 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is_ Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x ppp 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION : LAC IE Non-LACIE 

Trainin g y C lassification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE CLSNEW 

II. SOURCE (Include coder /person with knowledge of program usage): 

L.F. Guseman, Jr. / Bruce Marion 

III, PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Classify labeled measurement vectors using the B vector 

and associated decision regions previously determined by the 
minimization procedure 

V. ALGORITHM/ TECHNIQUE (short description; Include reference if available) 
Ref er enc e : LFSPMC 

Description: Each measurement vector is classified by first 
forming its inner product with the optimal B vector and then 
choosing the decision region of the result. An error matrix 
is built to record correct /incorrect classifications. 

VI. SPECIAL FEATURES: Each vector, its label, and assigned class may 
be output by setting appropriate option flags. 


VII, SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is not x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE ERRMAT 

II. SOURCE (include coder/person with knowledge of program usage) : 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Computes and outputs estimated confusion matrix based on the error 
matrix generated by the classification routine 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: LFSPMC 

Description: The ij th entry pij in the confusion matrix is the 
estimated probability that something in class j is classified 
into class i. 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x POP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Training Classificatio n x E stimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE RANCLN 

II. SOURCE (Include cader/person with knowledge of program usage): 

L.F. Guseman, Jr. / Bruce Marlon 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Generate and classify random vectors according to original 
class statistics to estimate probability of misclassif ication 
in the original feature space 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 
Reference: LFSPMC 

Description: Mean vector and covariance matrix for each class are 
used to produce random n-dimensional vectors. Each such vector 
is classified (Bayesian classifier) and a confusion matrix is 
estimated. 

VI. SPECIAL FEATURES: Variable dimensioned 


VII. SUBROUTINES REQUIRED: MCHLSK, GAUSS 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is_ Is no t x c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Trainin g x C lassificatio n y E stimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE MCHLSK 

II. SOURCE (include coder /person with knowledge of program usage): 
L.F. Guseman, Jr. / Bruce Marion 

I~I. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

1’V. PURPOSE: Compute the modified Cholesicy decomposition of the 
covariance matrix. 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: V.X. Van Rooy, M.S. Lynn, C.H. Snyder, ICSA Report No-275- 
025-008, Institute for Computer Services and Applications, 

Rice University, Houston, Texas . 

Description: Uses a modified method of Cholesky to decompose a 
positive definite symmetric matrix into upper and lower triangular 
matrices. 

VI. SPECIAL FEATURES: The decompositions overlay the elements of the 
original matrix. 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify . 

Software must be designed and develop for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x p pp 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE xNon-LACIE 

Trainin g x C lassificatio n x E stimation of Proportion 

x F eature Selection Signature Extension 

Temporal Samplin g S ystems Development Technique Development 

Crop Assessmen t , Y ield Estimation 

Other (specify) - 
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I. NAME: SUBROUTINE GAUSS 

II. SOURCE (include coder /person with knowledge of program usage): 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Compute a normally distributed random number of mean zero 
and covariance one. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: System/360 Scientific Subroutine Package, Version III, 

5th Edition Form GH20-0205-4, International Business Machines, 

August, 1970. 

Description: Twelve uniformly distributed random numbers on the 
interval 0 to 1 are added, giving a random number approximately 
normally distributed about 6; then 6 is subtracted, giving the 
desired result. 

VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: RANDU 

VIII. DOCUMENTATION: x A vailable Not Available Will be Available 

IX. _x_Software available off-the-shelf. Approximate Cost: 

Is Is not x compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non- LAC IE 

Training x Classification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) __ 
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I. NAME: SUBROUTINE RANDU 

II. SOURCE (include coder/person with knowledge of program usage) : 

L.F. Guseman, Jr. / Bruce Marion 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/65 

IV. PURPOSE: Compute a uniformly distributed random number between zero 
and one. 

V. ALGORITHM/TECHNIQUE (short description; Include reference if available) 

Reference: System/360 Scientific Subroutine Package, Version III, 

5th Edition, Form GH20-0205-4, International Business Machines, 

August, 1970. 

Description: 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

Is Is no t x compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Trainin g x C lassification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

C rop Assessment Yield Estimation 

Other (specify) 
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THIS FORM PREPARED BY LEC 

I. NAME: CANONICAL 

II. SOURCE (include coder /person with knowledge o£ program usage): 
Walter Eppler 

III. PROGRAMMING LANGUAGE : FORTRAN COMPUTER: IBM 360/75 

IV. PURPOSE: Too classify multispectral data. 


EPPLEH 
1 of 2 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See next page 

Description: See next page 


VI. SPECIAL FEATURES: See next page 


VII. SUBROUTINES REQUIRED: None 

VIII. DOCUMENTATION: Available Not Available x Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: Q 

I s xl s not compatible with Image 100 configuration in present form 

x Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

x T rainin g x C lassification Estimation of Proportion 

x F eature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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CANONICAL ANALYSIS FOR INCREASED CLASSIFICATION 
SPEED AND CHANNEL SELECTION 


Walter Eppler 

Lockheed Electronics Company 

Kenneth Baker 
NASA Johnson Space Center 


1 . Introduction 

t 

The well-known classification rule based on the maximum-likelihood criterion 
and assumed normal probability density functions involves evaluating quadratic 
forms in the case of M classes. Sometimes a linear 'transformation is performed 
on the original N measurements to form N < M measurements for use in evaluating 
the quadratic forms resulting in a reduction in computation time. The 
disadvantages of this approach are: 

1* Additional computer time is required to perform the linear transforma- 
tions . 

2. It is inevitable that some (usually small but unknown) class separa- 
bility is lost in the dimensionality-reduction. 

The algorithm described in this paper has the advantage that it uses only as 
many channels as necessary to make the desired discriminations; if necessary, 
all channels are used and no information is sacrificed. Classification speed 
is achieved by using first those channels which are most important for discrim- 
ination* The paper describes how the proper order is determined; the deriva- 
tion has obvious application in the general area of channel-selection. 

Empirical results are presented to show the optimum channel-order and resulting 
increase in classification speed in typical applications* 
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D. PATTERN RECOGNITION LIBRARY 

4. ESTIMATION OF PROPORTION PACKAGE 

A package which estimates proportions by using maximum likelihood, 
method of moments, estimated confusion matrix, and classification is 
being documented by Dr. Bill Coberly (currently at the University of 
Houston) . 
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I. NAME: SUBROUTINE APREST - ESTIMATION OF PROPORTION PACKAGE 

II. SOURCE (include coder /person with knowledge of program usage): 

L. F. Guseman, Jr. /Bruce Marion, Department of Mathematics, TAMU 

III. PROGRAMMING LANGUAGE : FORTRAN IV COMPUTER: 360/65 

IV. PURPOSE: To est ima te the proportion of each of m classes from a 
set of random observation vectors. 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: L. F. Guseman, Jr. and Bruce Marion, A method for estimating 
proportions, Contract NAS-9-13894-1S, Report #4, Dept, of Math., Texas 
A&M University, April 1975, 

Description: Given a classification procedure, its associated confusion 
matrix, and a set of random observation vectors selected from among m 
classes, the vector of proportions for the classes is estimated. 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: LSI, AWNW 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. _x Software available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpos e x P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Training Classificatio n x E st ima tion of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Y ield Estimation 

Other (specify) 


D.4-2 


COMPUTER PROGRAM IDENTIFICATION FORM 


2 of 9 


I. NAME; SUBROUTINE LSI 

II. SOURCE (Include coder/person with knowledge of program usage); 

L. F. Guseman, Jr. /Bruce P. Marlon, Dept. Math., TAMU 

III. PROGRAMMING LANGUAGE; FORTRAN IV COMPUTER; 360/65 

IV. PURPOSE: To solve the least squares problem with inequality constraints (LSI) 
minimize | |Ex-u| | subject to x > 0, £x^ = 1 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: c. L. Lawson and R. J. Hanson, Solving Least Squares Problems, 
Prentice-Hall, 1974, Chapter 23. 

Description: Converts LSI into least distance problem (LDP). 

VI. SPECIAL FEATURES: Checks for singular E (assumed square). 

VII. SUBROUTINES REQUIRED: SVDM, LDP 

VIII. D0C1P lATION: x Available Not Available Will be Available 

IX. jc Software available off-the-shelf. Approximate Cost; -?- 

I sx I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

G eneral purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

C rop Assessmen t Y ield Estimation 
Other (specify) 
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I. NAME: SUBROUTINE SVDM 

II. SOURCE (include coder /person with knowledge of program usage) : 

G. H. Golub, Computer Science Dept., Stanford Univ., Stanford, Ca. 94305 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 1108, 360/65, CDC6600 

ALGOL 60 

IV. PURPOSE: 

Singular value decomposition of rectangular matrix 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: g. H, Golub and C. Reinsch, "Singular value decomposition 

and least squares solutions," Numer. Math 14(1970), pp. 403-421. 

Description: 

Decompose a given matrix A into matrices U,S,V such that 
A = U * 3 * V T 

VI. SPECIAL FEATURES: 

Modification of SUBROUTINE DSVD. 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost; -7- 

I s x Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

_x General purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : x LACIEx Non-LACIE 

Trainin g x C lassificatio n x E stimation of Proportion 

x F eature Selectio n x Signature Extension 

Temporal Samplin g S ystems Development x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE LDP 

II. SOURCE (include coder /person with knowledge of program usage): 

L. F. Guseman, Jr. /Bruce D. Marion, Dept, of Math,, TAMU 

III. PROGRAMMING LANGUAGE : FORTRAN IV COMPUTER:360/65 

IV, PURPOSE: To solve the least distance problem (LDP): minimize |jg|| 
subject to Ag > w 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See LSI 

Description: Converts LDP into nonnegative least squares problem 
(NNLS) . 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: NNLS 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost: 

I s x Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than i-ORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

y G eneral purpose Special purpos e P rogram package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non- LAC IE 

Training Classificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

O ther (specify) 
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I. NAME: SUBROUTINE NNLS 

II. SOURCE (Include coder/person with knowledge of program usage): 

L. F. Guseman, Jr. /Bruce P. Marion, Dept. Math., TAMU 

III. PROGRAMMING LANGUAGE : FORTRAN IV COMPUTER: 360/65 

IV. PURPOSE: To solve the nonnegative least squares problem (NNLS): 
minimize | | Qs-v j | subject to z > 0 . 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See LSI 

Description: Given an m><n matrix C, and an m-vector v, compute an 
n-vector z with nonnegative entries which minimizes 
| j Cz-v 1 1 . 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: H12, Gl, G2, DIFF 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX‘. _x Software available off-the-shelf. Approximate Cost: 

l ax I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify ___ 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

x G eneral purpose Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : y LAC IE x Non-LACIE 

Training Classificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) _ ________ 
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I. NAME: SUBROUTINE H12 

II. SOURCE (include coder/person with knowledge of program usage): 

L. F. Guseman, Jr. /Bruce P. Marion, Dept. Math., TAMU 

III. PROGRAMMING LANGUAGE : FORTRAN IV COMPUTER: 360/65 

IV. PURPOSE: construct and apply a Householder Transformation: 

Q = I + U(U T )/B 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference : see LSI 


Description* Given m-vector v and integers 1^ and 1^, computes 

m-vector u and S such that symmetric orthogonal matrix 
A = I + (uu 2 ’)/(SU 1 ) satisfies Q v - w. 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

__ In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION :x LACIEx Non-LACIE 

Training Classificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development_x_Technique Development 

Crop Assessment Yield Estimation 

Other (specify) _____ 
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I. NAME: SUBROUTINE G1 

II. SOURCE (include coder/person with knowledge of program usage): 

L, F. Guseman, Jr. /Bruce P. Marion, Dept, of Math., TAMU 

III. PROGRAMMING LANGUAGE : FORTRAN IV COMPUTER: 360/6? 

IV. PURPOSE: Construction and application of rotation matrices 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: see LSI 

Description: Given c, s, z^, computes 

■(-: :) (3 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LAC IE x Non-LACIE 

Training Classificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE G?. 

II. SOURCE (include coder /person with knowledge of program usage): 

L. F. Guseman, Jr. /Bruce f. Marion, Dept. Math., TAMU 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 350 /65 

IV. PURPOSE: Construction and application of rotation matrices. 

V. ALGORITHM/TECHNIO' T E (short description; include reference if available) 
Reference: See LSI 

Description: Applies the rotation computed by G1 to x^, x 2 : 



VI. SPECIAL FEATURES: 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. _x Software available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI . APPLICATION: x LACIEx Non-LACIE 

Training Classificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: SUBROUTINE DIFF 

II. SOURCE (include coder/person with knowledge of program usage) : 

L. F. Guseman, Jr. /Bruce P. Marion, Dept. Math., TAMU 

III. PROGRAMMING LANGUAGE : FORTRAN C0MPUTER:360/65 

IV. PURPOSE: To test n IF((x+h)-x) f 0 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: See LSI 

Description: The intermediate sum 2 = x+h is computed in the calling 

program using n precision arithmetic. The difference 2 -x 
is computed by this function using n precision arithmetic. 

VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x Software available off-the-shelf. Approximate Cost: 

I s x I s not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

x PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: x LACIE x Non-LACIE 

Training Classificatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t xT echnique Development 

Crop Assessment Yield Estimation 

. Other (specify) 


jT 
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I. NAME: PROPORTION ESTIMATION PACKAGE 

II. SOURCE (include coder/person with knowledge of program usage): 

W. A. Coberly, Dept, of Math. Sci., U. of Tulsa 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER:360/75 

IV. PURPOSE: Compute proportion estimates for use in acreage estimation. 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: "An empirical comparison of five proportion estimates," 

Coberly and Odell, NASA/ JSC Tech. Report, 1974, U. of 
Texas at Dallas. 

Description! 


VI. SPECIAL FEATURES: Five seperate driver programs 


VII. 


VIII. 


St. JTINES REQUIRED: 

*CLASS *MIXTUR 
Main *MXLIKE *MOMENT 
*CONF 

DOCUMENTATION: X Available 


MINV 

Subroutines PREP 

QUADPR 

Not Available Will be Available 


IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 


X. IMPLEMENTATION LEVEL 

General purpos e X Special purpos e X P rogram package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: LACIE Non-LACIE 

T rainin g C lassificatio n X E stimation of Proportion 
Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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DATA PRODUCTS REQUIRED 

1. Unlabeled (recognition) data set in intermediate format (vector by 
vector) . 

2. Test data set in above format, 

3. Labeled (training) statistics, 

4. Histograms of labeled data. 

5. Unlabeled statistics. 

6* Histograms of unlabeled data. 
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DATA FLOW 

1. Classification estimator (a^) • 



3. Odell-ChHikara confusion matrix estimate (a^) 
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4. Minimum Chi-square estimate (ct^.) 


5. Moment estimate (c^jq) 
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I. NAME: CLASS 

II. SOURCE (Include coder/ person with knowledge of program usage): 

W. A. Coberly, Dept, of Math., U. of Tulsa 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: 360/75 

IV. PURPOSE: Classify observation vectors (pixels) 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: "An empirical comparison of five proportion estimates," 

Coberly and Odell, NASA/JSC Tech. Report 1974, U. of 
Texas at Dallas 

Description: 

Maximum likelihood classification rule 

VI. SPECIAL FEATURES: Inputs are training signatures and unlabeled data. 

The likelihoods may be saved on temporary disk 
storage. 

VII. SUBROUTINES REQUIRED: MINV 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. x S oftware available off-the-shelf. Approximate Cost; 

Is Is no t c ompatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e X Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : X LAC IE X Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

O ther (specify) 
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I. NAME: MXLIKE 

II. SOURCE (include coder/person with knowledge of program usage): 

W. A. Coberly, Dept, of Math., U. of Tulsa 

III. PROGRAMMING LANGUAGE; FORTRAN IV COMPUTER: IBM 360/75 

IV. PURPOSE: Maximum likelihood solution of proportions. (Fixed 

signatures) 

V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 

Reference: "The Numerical evaluations of maximum-likelihood 

estimate of mixture proportions," Peters and Coberly 
NASA/ JSC Tech. Report, 1975, UTD. 

Description: 

Iteration scheme for finding maximum likelihood 
estimates of proportions, 

VI. SPECIAL FEATURES: Uses output of likelihoods for each 

from CLASS, stored in a temporary disk file. 


VII. SUBROUTINES REQUIRED: NONE 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

' X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION : X LAC IE Non-LACIE 

Training Classification Estimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t X T echnique Development 

C rop Assessment Yield Estimation 

Other (specify) _____________________________________ 
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I. NAME: CONF 

II. SOURCE (include coder /person with knowledge of program usage): 

W. A. Coberly, Dept, of Math., U. of Tulsa 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/75 

IV. PURPOSE: Proportion estimation by Odell-Chhikara method using the 

confusion matrix. 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: "An empirical comparison of five proportion estimates," 

Coberly and Odell, NASA/ JSC Tech. Report, 1974, UTD. 

Description: Finds constrained solution of the equation 

a = Pa where & is the classification estimate of 
the proportions and P is the confusion matrix. 

VI. SPECIAL FEATURES: The confusion matrix P is estimated by the 

classification of test data. 


VII. SUBROUTINES REQUIRED: QUADPR, PREP 

VIII. DOCUMENTATION : X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e X Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: * LACIE Non-LACIE 

Training Classi r lcatio n x E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development x T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. NAME: MIXTUR 

II, SOURCE (include coder/person with knowledge of program usage): 

W. A. Coberly, Dept, of Math., U. of Tulsa 

III. PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/75 

IV. PURPOSE: Proportion estimation by method of minimum chi-square. 


V. ALGORITHM/ TECHNIQUE (short description; include reference if available) 

Reference: "An empirical comparison of five proportion estimates," 

Coberly, Odell, NASA/JSC, Tech. Report, 1974, UTD . 

Description: Constrained least squares solution of system found by 
equating the percentiles of the marginal mixture 
distributions . 


VI. SPECIAL FEATURES: 


VII. SUBROUTINES REQUIRED: QUADPR, PREP 

VIII. DOCUMENTATION: X A vailable Not Available Will be Available 

IX. X Software available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

General purpos e x Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION :__^LACIE Non-LACIE 

Training Classificatio n X E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Developmen t ^ T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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II. 

III. 

IV. 

V. 


VI. 

VII. 

III. 

IX. 


X. 

XI. 
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NAME: M0MENT 

SOURCE (Include coder/ person with knowledge of program usage) : 
U. A. Coberly, Math. Dept, U. of Tulsa 

PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IBM 360/75 

PURPOSE: Proportion estimation by method of moments. 


ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: "An empirical comparison of five proportion estimates," 
Coberly, Odell, NASA Tech. Report 1974, U. of Texas at 
Dallas 

Description: Constrained weighted least square solution of mixture 

moment equations. 


SPECIAL FEATURES: 


SUBROUTINES REQUIRED: QUADPR, PREP 

DOCUMENTATION: X Available Not Available Will be Available 

Software available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

IMPLEMENTATION LEVEL 

General purpos e ^ Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

APPLICATION: X LACIE Non-LACIE 

Training Classificatio n X E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development X T echnique Development 

Crop Assessment Y ield Estimation 

Other (specify) 
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I. NAMES QUADPR 

II. SOURCE (include coder/person with knowledge of program usage): 
W. A. Coberly, Dept, of Math., U. of Tulsa 

III. PROGRAMMING LANGUAGE: FORTRAN I? COMPUTER: IBM 360/75 

IV. PURPOSE: Constrained quadratic programming scheme (Wolfe) 


V. ALGORITHM/TECHNIQUE (short description; include reference if available) 

Reference: "Optimization techniques with FORTRAN," Kuester and 

Mize, McGraw Hill, pp. 106-119. 

Description: Quadratic programming routine used to solve constrained 
least squares problems. 


VI. SPECIAL FEATURES: Does not handle singular quadratic objective 

functions. 


VII. SUBROUTINES REQUIRED: None 

VIII. DOCUMENTATION: X Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

X G eneral purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: X LAC IE Non-LACIE 

Training Classificatio n X E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development T echnique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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I. 

II. 

III. 

IV. 

V. 


VI. 

VII. 

VIII. 

IX. 


X. 

XI. 
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NAME: PREP 

SOURCE (include coder/person with knowledge of program usage) : 

W, A. Coberly, Dept, of Math., U. of Tulsa 

PROGRAMMING LANGUAGE: FORTRAN IV COMPUTER: IB M 360/75 

PURPOSE: Prepare objective function from a system of linear equations 

for input to QUADPR. 


ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: 


Description: 


Constructs quadratic objective function for constrained 
least square solution by Wolfe’s method. 


SPECIAL FEATURES: 


SUBROUTINES REQUIRED: NONE 

DOCUMENTATION: x A vailabl e N ot Available Will be Available 

Software available off-the-shelf. Approximate Cost: 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

I n language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

IMPLEMENTATION LEVEL 

X G eneral purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

APPLICATION: x LACIE Non-LACIE 

T rainin g C lassificatio n x E stimation of Proportion 
Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

.Other (specify) 
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I. NAME: MINV 

II. SOURCE (Include coder /person with knowledge of program usage): 

W. A. Coberly, Dept, of Math., U. of Tulsa 

III. PROGRAMMING LANGUAGE: FORTRAN COMPUTER: 360/75 

IV. PURPOSE: Invert a matrix 

V. ALGORITHM/TECHNIQUE (short description; include reference if available) 
Reference: IM SS p_xxx Manual 

Description* An elementary Gaussian elimination routine 

VI. SPECIAL FEATURES: a. routine inverting only a symmetric positive- 

definite covariance matrix is required. 

VII. SUBROUTINES REQUIRED: 

VIII. DOCUMENTATION: x Available Not Available Will be Available 

IX. X S oftware available off-the-shelf. Approximate Cost: -0- 

Is Is not compatible with Image 100 configuration in present form 

Non-PDP 11/45 FORTRAN 

Some machine language conversion necessary 

In language other than FORTRAN. Specify 

Software must be designed and developed for Image 100 

X. IMPLEMENTATION LEVEL 

X G eneral purpose Special purpose Program package 

PDP 11/45 FORTRAN LIBRARY FUNCTION OR SUBROUTINE 

XI. APPLICATION: XLACIE Non-LACIE 

Training Classificatio n X E stimation of Proportion 

Feature Selection Signature Extension 

Temporal Sampling Systems Development Technique Development 

Crop Assessment Yield Estimation 

Other (specify) 
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D. PATTERN RECOGNITION LIBRARY 

5. SPATIAL PATTERN RECOGNITION PACKAGE 


D. 5-1 


D. PATTERN RECOGNITION LIBRARY 


6. SIGNATURE EXTENSION PACKAGE 

A program is being developed at the University of Houston to 
perform signature extension using the maximum likelihood function. A 
similar, but more specialized, algorithm for signature extension is 
being implemented by LEC. A specialized version of the LEC algorithm 
is being studied by ERIM (MASC Algorithm). An additional signature 
extension technique is being developed by Dr. Ungar at GFSC (no details 
available) . 
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IMPLEMENTATION PRIORITY SUMMARY 


The following is a summary by package of those subroutines and 
programs for which identification forms have been completed. Each 
program was assigned one of the implementation priorities defined in 
Section IV and repeated below for ease of reference. 

Priority I - Basic utility subroutine or program 
required for the implementation of 
existing image analysis techniques 
Priority II - Existing image analysis techniques 

(requires majority of Priority I items) 

Priority III - Subroutine or program required to 

support development and/or refinement 
of image analysis techniques 

Priority IV - Image analysis techniques in research and 
development stage. 
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PROGRAM 


PRIORITY 


A. PRE-ANALYSIS DATA PROCESSING LIBRARY 

1. PREPROCESSING PACKAGE 


GENERAL RADIOMETRIC CORRECTION (G.E.) I 

2. REGISTRATION PACKAGE 

SUBROUTINE CORLAT II 

SUBROUTINE REGSTR II 

3. TRAINING FIELD SELECTION PACKAGE 

IRREGULARLY-SHAPED FIELDS I 

LINE FIND ROUTINE I 

GREY SCALE MAP PACKAGE IV 

KONTUR PACKAGE IV 


4. SIMULATED DATA GENERATOR PACKAGE 


PROGRAM PRIORITY 

B. GENERAL MATHEMATICAL LIBRARY 

1. LINEAR ALGEBRA PACKAGE 

SUBROUTINE MINVSP I 

SUBROUTINE MINVDP I 

LU-Decomposltion Package I 

SUBROUTINE GINV2M I 

SUBROUTINE CHLSICY I 

TSTMIS (single precision) I 

TSTMID (double precision) I 

SUBROUTINE SEIGEN I 

SUBROUTINE DSVD I 

SUBROUTINE MATMUL I 

SUBROUTINE EIGEN I 

2. OPTIMIZATION PACKAGE 

SUBROUTINE DFMFP I 

SUBROUTINE FRKWLF I 

SUBROUTINE QM431 I 

SUBROUTINE LSI I 

3. APPROXIMATION PACKAGE 

Basic Spline Routines — Single Variable I 

4. TRANSFORM PACKAGE 

(MSFAR) Multispectral Fourier Analysis Routines III 
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PROGRAM 

C. GENERAL STATISTICAL LIBRARY 

1. STATISTICAL SUMMARY PACKAGE 

SUBROUTINE MECOV 
SUBROUTINE GROP1 
SUBROUTINE MEDIAN 
SUBROUTINE MODE 
SUBROUTINE RANGE 
SUBROUTINE MNDEV 
SUBROUTINE SDEV 
SUBROUTINE ORDER 
PLOT PAC 

BINARY SEARCH/N-DIMENSLON HISTOGRAM ACQUISITION (G.E.) 

2. DENSITIES & DISTRIBUTIONS PACKAGE 

FUNCTION PHINV (P) 

DOUBLE PRECISION FUNCTION DPHI 
REAL FUNCTION FISH 
DOUBLE PRECISION FUNCTION DERF 
SINGLE PRECISION FUNCTION ERF 
DOUBLE PRECISION FUNCTION DERFC 
SINGLE PRECISION FUNCTION ERFC 
REAL FUNCTION CHI 2 
FUNCTION EXPAN 
REAL FUNCTION POSI 
FUNCTION FISHIN 
REAL FUNCTION BIN 
FUNCTION HYTRIC 
REAL FUNCTION BET INC 
REAL FUNCTION GAMIN 
REAL FUNCTION STUD IN 
REAL FUNCTION ZOT 
FUNCTION RANDOM 
REAL FUNCTION ZIP 
REAL FUNCTION ZARF 
FUNCTION KPSNSM 
SUBROUTINE BIN 2 
SUBROUTINE OVERFL(I) 

3. REGRESSION PACKAGE 

4. STATISTICAL TESTS 

SUBROUTINE TTEST 
SUBROUTINE C0VARI 
SUBROUTINE ANOV1 
SUBROUTINE ANOVA2 
SUBROUTINE RBANOV 
FUNCTION ZOR 
SUBROUTINE KOLSMR 
SUBROUTINE FACTOR 


PRIORITY 


I 

I 

I 

I 

I 

I 

I 

I 

I 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 
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I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


I 

I 

I 

I 

I 

I 

I 

I 
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PROGRAM 


PRIORITY 


D. PATTERN RECOGNITION LIBRARY 

1. CLUSTERING PACKAGE 

SEMI-SUPERVISED CLUSTER ANALYSIS (G.E.) II 

ISOCLS Routines II 

2. CLASSIFICATION PACKAGE 

ELLTAB I 

MAXIMUM LIKELIHOOD RESOLUTION (G.E.) I 

3. FEATURE SELECTION PACKAGE 

UH Linear Feature Selection II 

LF MC : Linear Feature Selection (PMC) II 

CANONICAL II 

4. ESTIMATION OF PROPORTION PACKAGE 

SUBROUTINE APREST II 

5. SPATIAL PATTERN RECOGNITION PACKAGE 


6. SIGNATURE EXTENSION PACKAGE 


